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ELECTROMAGNETIC INDUCTION IN A SEMI-INFINITE 
CONDUCTOR WITH A PLANE BOUNDARY 


By A. T. PRICE (Imperial College, London) 


keceived 14 March 1950] 


SUMMARY 
The general theory of the induction and of the free decay of electric currents in 
a semi-infinite conductor with a plane boundary is considered. The solutions of the 
fundamental equations are of two types, which correspond to two distinct modes 
of decay of currents in the conductor. The properties of these modes are investigated. 


Detailed calculations are made of the induced field and currents produced by a 


periodic or ape iodic line current parallel to the face of the conductor, and the method 
of calculation for any other inducing fields is indicated. An explanation is given of 
the indeterminate character of certain problems. 


1. Introduction 

Tue theory of the electric currents and magnetization induced in a 
permeable conducting sphere by a varying magnetic field has been the 
subject of numerous investigations by Lamb (1), Price (2), Chapman and 
Price (3), Lahiri and Price (4), Elsasser (5), and in many earlier references 
quoted in (4). The main object of these investigations has been to obtain 
results applicable to various geophysical problems (6). These results are 
expressed in terms of spherical harmonics, so that a necessary condition 
for the practical application of the theory is that the field considered shall 
be expressible by a series of spherical harmonics which is fairly rapidly 
convergent. This rules out the possibility of considering in this way fields 
of decidedly local character, since these would involve many harmonies of 
high order. It is obviously possible, however, to obtain some information 
in such cases by treating the earth as a semi-infinite conductor bounded 
by a plane. 

A number of investigations relating to such a semi-infinite conductor 
have been made by Niven (7), Sommerfeld (8), Slichter (9), Buchholz (10), 
but these have usually dealt only with some specific problem, or have been 
specially concerned with the radiation field; a general theory, similar to 
that developed for the spherical conductor and dealing specially with the 
induction field (as distinct from the radiation field treated by Sommerfeld 
and others), does not seem to have been fully considered. This theory is 
developed in the present paper and the physical interpretation of the two 
distinct types of fundamental solutions which are found of the relevant 
equations is examined. The methods of calculating the induced field and 
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currents corresponding to any inducing magnetic field of external origin 
are also indicated. 

The particular case when the inducing field is produced by a periodic 
current along a line parallel to the face of the conductor is worked out in 
detail, and numerical results are obtained which permit the induced field 
to be found easily for any line current situated at any distance from a 
conductor of uniform conductivity. The cases of a circular current circuit 
and of an oscillating magnetic dipole are dealt with by A. N. Gordon in 
another paper (11). 

In connexion with geophysical applications of the theory, it is of course 
obvious that it is not always possible to infer the induced field over a small 
region on the surface of a conductor (e.g. a restricted region of the earth’s 
surface) from a knowledge of the inducing field over only the same limited 
region. This is further emphasized in the present theory by the fact that 
certain induction problems relating to the semi-infinite conductor are 
indeterminate. These are considered in section 12 and the corresponding 
restriction on the applicability of the theory to geophysical problems is 
indicated. 

When the inducing field is not periodic, any change in its intensity will 
induce a distribution of.currents in the conductor which subsequently 
diffuse through the conductor and decay. In the case of a sudden change 
of the field, the initial distribution of currents will be confined to a thin 
layer near the surface, but will gradually penetrate to greater depths. 
This case is worked out in sections 13-16; the result for any other time 
variation of the field can be obtained from this particular case by a well- 
known formula. 

Since it is necessary to make frequent references to the fundamental 
field equations and boundary conditions, these are briefly summarized in 
sections 3 and 4. 


2. Notation 


The following notation is used throughout: 


H, magnetic field intensity (, magnetic scalar potential 
B, magnetic induction ¢, electric potential 

i, conduction current density p, volume density of charge 
E, electric field intensity o, surface density of charge 
D, electric displacement u, Magnetic permeability 

I, intensity of magnetization Kk, electric conductivity 

A, magnetic vector potential e, dielectric constant 


All quantities are measured in Gaussian units. 
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igin 3. General equations of the field 
Inside a continuous medium the field vectors satisfy Maxwell’s equations 
clic l 
t in curl H (47i+-D), (3.1) 
} r 
field 
ma curl E B, (3.2) 
° Cc 
cult 
i togethe1 with 
n in —_ : 
divB = 0 div i p, div D = 4p. (3.3) 
urse If the medium is isotropic and uniform, we have also 
mall | B = pH i = «Bb, D = «€E, (3.4) 
rth’s where », «, and € are constants. In this case it follows from (3.1) that 
ited | 
that div(47«,«)D+div D 0, 
* are whence p Po ' Ake) t (3.5) 
\ding fos pm ee ee ' 
showing that any initial volume charge distribution in the conductor dis- 
ms 1s . . 4s . 
perses at a rate uninfluenced by other field conditions. Hence there is no 
a oss of generality in taking p zero, and consequently 
y will 
iv E = 0. 3.6 
ently div I —) 
lange It then follows, on taking the curl of equations (3.1) and (3.2), that all the 
thin omponents of every field vector satisfy an equation of the form 
»pths. — UU c e2 : on 
as V7 4 —| 47K + €--~ F. (3.7) 
time Cc ct ot*| 
well , : ‘ 
In a conductor, when the time changes in the field are not too rapid, 
the term ecF ct may be ignored compared with 4z7«F. The equation 
enta ae ; . ogee ‘ 
Te 3.7) then reduces to the induction (or diffusion) equation 
zed in 
eae 47K OC 7, —— 
V7F . F. (3.8) 
co a 
In a dielectric « is zero or negligibly small, so that (3.7) reduces to the 
ave equation 2 
= te C7 F , 
VF =O, (3.9) 
ial CC oF 
howing that electromagnetic effects are propagated with speed c/,/(eu) 
rge : through the medium. When we are concerned only with the induction field 
rge rt relatively slow field-variations, the time of travel across any portion 


| ofthe dielectric, within the region considered, is negligibly small compared 

ith the time taken by the field changes at a point to become effective. 

Hence we can ignore the term on the right of (3.9), so that in the dielectric 
e have simply 


V2F = 0. (3.10) 
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The above approximations are, of course, equivalent to ignoring the 
magnetic effects of the displacement current both in the conductor and in 
the dielectric. It follows from (3.1) that, in the dielectric, curl H is zero 
to the same order of approximation. Hence H can be derived (approxi- 
mately) from a scalar potential, say 


H grad Q), (3.11) 


where © satisfies Laplace's equation since div H is zero. 
If H and E are derived from the vector and scalar potentials (A, 4) of 
the field by the usual relations, 
wH curl A, E A c—grad ¢ (3.12) 
it is found that both A and ¢ satisfy the fundamental equation (3.7) (or 
the corresponding approximate equations (3.8) or (3.10)) provided A and 
¢ satisfy the relation 
° AL c 
div A! (42 1 € - \¢ 0. (3.13) 
"r ct 
The results obtained in this paper were originally expressed in terms of 
A and 4, but Dr. A. N. Gordon points out that it is somewhat simpler to 
obtain them directly in terms of E and H. This is because E and H are 


non-divergent, whereas A in general is not. 


4. Boundary conditions 

At the interface between two media, the tangential components of E and 
of H, and the normal component of B, are continuous. The corresponding 
boundary conditions for A and ¢ are that the tangential components of A 
and of grad¢ are continuous. The normal components of D and of i are 
not in general continuous, the sum of their components in the outward 
directions from the interface being 470 and —6é, respectively. 

The determination of the currents induced in a conductor of any shape 
by a given external field, or of the free decay of a system of currents in 
the conductor, involves finding appropriate solutions of equations (3.8) and 
(3.10) for the components of the field vectors, which must also satisfy the 


above boundary conditions. 


5. Semi-infinite conductor with a plane boundary 
Taking rectangular axes, suppose the conductor occupies the half-space 
z <0, so that its boundary is the plane z = 0, and suppose a varying 


magnetic system, consisting of moving magnetic matter or of electric 


current systems, is situated in the region z > h > 0. In the intermediate 
region 0 2 h we assume that the medium is non-conducting and 


non-magnetic. 
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The elementary solutions of the differential equations and boundary 
mditions of sections 3 and 4 will be shown to be of two types. In the 
first type the corresponding current systems in the conductor possess a 
magnetic field outside it: in the second they do not. When we are con- 
erned solely with the effects produced by a varying external magnetic 
system, we require solutions of the first type only. For it will be shown 
that a solution satisfying the differential equations and boundary con- 
litions and containing a part representing the inducing field can be 
obtained in terms of solutions of the first type. Now if another solution 
is found also consistent with the inducing field, the difference between 
the two solutions would also be a solution (since all the equations are 
linear) and moreover would correspond to a zero inducing field. Hence 
this difference can only represent a freely decaying system of currents in 
the conductor. It would consist, in general, of solutions of both types, but 
vould necessarily be zero if the initial current intensity were everywhere 


reToO 


6. The elementary solutions of the equations 

Since the tangential components of E and H are continuous at the 
boundary z — 0, while the conductivity « and permeability are discon 
tinuous there, it is convenient to express each of these field vectors in 
terms of elementary Solutions of the form Z(z,/)F(a,y), the variable ¢ 
ecurring only in Z because @ ct is multiplied by «xp in the relevant 


equation (3.8) 


Thus, taking E Z(z, t) F(x, y) (6.1) 
ind substituting in (3.6), (3.8), and (3.10), we obtain 
7 (OF. F ,OZ 
Z\° ; “| F‘ 0, (6.2) 
Ca cy " OZ 
ceF Cc? l [4a Z ; =e 
. ae C »)E tor Z 0 (6 3) 
Ca cy- Z » ct Cz* 
c’eF oF lo. -« 
ne . . as —-S for z> 0. (6.4) 
os" cy” Z 02? 
Since the coefficient of F in the term on the right of (6.3) or of (6.4) is 
ndependent of « and y, it must be a real constant,* —A* say, so that 
CF oF _ 
-+. 2 0). (6.5) 
Ca cy~ 


it this constant is necessarily negative. 
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eZ a dak CZ 
and — — A?Z ct 0 for 2z (0). (6.6) 
02° c* oat 
PE iss ' 
-—A°Z=0 for z>0. (6.7) 
22 
From (6.2) it follows that either 
, F. F 
F. QO and cd pci 0. (6.8) 
d CX cy 
1 (oF, oF, 1 cZ d 
or a + "| - X (6.9 
Fi\ cx oy Z oz 


where « is any constant (real or complex). 
The two cases corresponding to (6.8) and (6.9), respectively. will now 


be considered separately. 


6.1. Elementary solutions of the first type 


In the first case, it is easily seen from (6.8) that F can be written in the 


form E éP j P 


0). (6.10) 
cy Cr 


where P is a function of x and y, which from (6.5) must satisfy the equation 
OP 8 5 “s 

a ~+A*P 0. (6.11) 

ox cy~ 


In the dielectric we have. from (6.7) 


Z A (t)e** 4 BY (t)je-“=, 


. ; . ; \=)( ¢ rP oP : 
so that E iJ (the“+t Bite ey , 0) for 2 0, (6.12) 
cy ox 
and, since the tangential components of E are continuous at 2 0. the 


corresponding expression for E inside the conductor is 
€ P C P 


cy Cr 


E Z(z,t)| 0] for 2z 0. (6.13) 


where Z(—0, t) A(t) +-Bt). (6.14) 


The corresponding magnetic field is readily obtained from equation (3.2), 
which shows that 
: cZceP cZeP 
wH c ; 


OZ OX oO2 OY 


,A°ZP for 2 0. (6.15) 


and H cA grad{(.A/(t)e=— A(the-)P! for z 0. (6,16) 


These expressions satisfy the boundary condition that the normal com- 
ponent of B (or B) is continuous at z = 0 in virtue of equation (6.14) 
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The remaining boundary condition, that the tangential components of H 
are continuous, leads to the further relation 

Zr {—®; 2) Au}.A(t)—B(t)}. (6.17) 
where Z’ denotes ¢Z/éz. 

The above expressions for the field vectors will be called solutions of the 
first type, and it will now be shown that only solutions of this type are 
required when considering the effects of an external inducing magnetic 
field. For in this case the field intensity in the conductor must tend to 
ero as z becomes large and negative. so that 

Z(z,t)>0 as 2>—0o. (6.18) 

Further, the scalar potential Q of the total magnetic field will satisfy 

Laplace’s equation in the region 0 Z h, and can therefore be expressed 
is @ SuMMation of functions of the form 

Q) YA (t)e* B(tye A=! P(x, y), (6.19) 

here A is real and positive (or zero*) and P satisfies equation (6.11). 

The function P(x, y) can, in fact, be written in polar coordinates in the 
form J.(AR)cos(s@+-e), and an arbitrary function 2, given over z = 0, can 
ye expressed in terms of these functions by means of Fourier’s theorem 
ind the Fourier—Bessel integral. For example, the potential at (x,y,z) of 


magnetic pole at (29, Yo, Z)) can be expressed in the region 0 < z < 2, in 
the form —m | e*~J,(AR) dA, where R = .{(x—2x9)?+(y—Yo)*}- 

Now the part of 2 in (6.19) which involves e* corresponds to sources 
f the field over the region z > h, i.e. to the given inducing field, while 
the part involving e~““ corresponds to the field arising in the region z < 0, 

to the field of the induced currents (and induced magnetism if u + 1) 


the conductor. Hence we may regard A(f) in (6.19) as a known function, 

nd on comparing with (6.16) we see that 
cA.oA (t) {(t), cAB(t) = Bit). (6.20) 
that .o7(t) is known. 

The equations (6.14), (6.17), and (6.18), together with (6.6) then serve 
to determine A(t) and Z(z,t) uniquely in terms of the given A(t). The 
nduced field will therefore be completely determined if the initial con- 
lition (i.e. the initial value of B(t)) is known. The induced field (and 
urent distribution) is thus expressed entirely in terms of solutions of the 
tst type. In these solutions the z-component of E, and therefore also of i, 
s everywhere zero. We have therefore the important result that th 
duced currents flow everywhere parallel to the surface of the conductor. 


* The case 0 requires special treatment and is considered in section 11. 
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7. Elementary solutions of the second type 
The second of the two possibilities mentioned in section 6, namely, the 
case when equations (6.9) are satisfied, will now be considered. These 
equations show that Z must be of the form 
Z C(t)e™. (7.1) 


Hence. inside the conductor we have. from (6.6), 


Sark ‘ 
(a2—A2)C(t) “EC(t), 
72 
so that C(t) = Cye-Pr, (7.2) 
2 A2 2 
where B = : ) (7.3) 
dirk 


From (6.9) and (6.5) it is easily shown that 

F Mie F, ; OP, F vx Fr, OP, 
~ 2 Ox” Oy [ oy 
where P, is a function of x and y satisfying (6.11). The terms involving P, 
obviously lead again to solutions of the first type and can therefore be 
omitted. Also from (6.5) the component F, satisfies (6.11). Hence, writing 


F. = AP, we have a second type of solution for E of the form 
(.. oP 1p 
E = C,e™ (2S OF yp) for 2 <0. (7.4) 
A cx Acy 


Outside the conductor, Z satisfies the equation (6.7) and on comparing 
this with (7.1) we see that « _r for z > 0. It then follows that E is 
of the form 


E = grad{A,(t)e*-4 B,(te“}P(x,y) for z> 0. (7.5) 
For the corresponding magnetic field inside the conductor we have 
. oP oP , 
wH ceurlE eC, e~Bl+02(y 2) , ——, 0} for 2z< 0. 
cy ox 
(7.6) 


Outside the conductor, since E is a gradient vector, we have 


H=0 for z>0. (7.7) 


Hence, the varying current distribution associated with the second type of 


solution has zero magnetic field outside the conductor. This explains why 
this type of solution does not enter into the problem of finding the currents 
induced by a varying external magnetic field. 

With regard to the boundary conditions, we first note that the con- 
tinuity of the normal component of B is already ensured since this 
component is everywhere zero. To make the tangential components of H 
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ontinuous (i.e. zero at z = 0) we should require C, to be zero if « were real. 
Hence, to get a non-trivial result we must take a and (, complex, and then 
make the 
A,+iB, and a d+iw, it is easily found that this condition 


eal part of C,e Pe %(A—a?/A) zero at z= 0. On writing 


0 
requires Ay = 0 and d = 0. It then follows that f is real and positive, 
nd, on taking the real parts of the expressions (7.4) and (7.6) for E and 
H. we obtain 

w cP w eP 


E By COS We COS Wz AP sin wz] for o* 0), (7.3) 
A Ox A cy 
' cP oP " 
uH cB,(A+w* Aye SIN w2 P : 0) for 2z 0. (7.9) 
cy Cx 
hese clearly satisfy div E — 0, divH = 0, which are equivalent to (6.9). 


One further boundary condition remains to be satisfied, namely, that 
the tangential components of E are continuous. Using (7.5) and (7.8), this 


gives the relation a 
A,(t)+ B,(t) Bye Pad. (7.10) 


It will be observed that A,(f) corresponds to an electric field of external 
rigin, and B,(¢t) to the field of a charge distribution on the surface of the 
onductor, since there is no volume charge in virtue of (3.5) and (3.6). Now 
if we suppose the initial conditions (at ¢ = 0, say) are given, all the 
parameters appearing in the expressions (7.8) and (7.9) will be determined, 
independently of A,(f). In particular the time factor e-f will be indepen- 
lent of A,(f), the equation (7.10) merely serving to determine the surface 
charge corresponding to B,(f). It follows that the field and current distri- 
bution within the conductor decay exponentially without change of form, 
ind are una ffe cted by any external ele ctric fie ld. 

The reason why this result holds is that we have assumed at the outset 
that the magnetic effects of the displacement current D can be ignored. 
hese effects are only important when the field is changing sufficiently 
rapidly to make cE comparable with 47xE, and to find them it would, of 
ourse, be necessary to use the general equation (3.1). 

The solutions of the second type which have been obtained above 
bviously correspond to the decay of current-distributions of special form 
in the conductor. The decay of current distributions in general will now 


e considered 


8. The free modes of decay of current in the conductor 

The results of the last section show that solutions of the second type 
ilwavs correspond to the free decay of certain distributions of electric 
‘urrent in the conductor, these distributions remaining unchanged in form 


vhile decreasing exponentially in intensity. There are also freely decaying 
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current systems corresponding to the first type of solution, which similarly 
remain unchanged in form while decreasing in intensity. These are ob 
tained by putting A(f) (corresponding to external magnetic sources) zero 
in (6.19). This implies that .o7(t) given by (6.20) is also zero; the boundary 
conditions (6.14), (6.17) then lead to 
Z'(0, t) AuZ (0, t). (8.1) 
Also the solutions of (6.6) which are of the required form f(z)g(t) are 
easily found to be 


z; (a, ¢ uz i b¢ x2 Ye bt (8.2 


where f is again given by (7.3), and «a is a constant, which, in order to 
make the time factor g(f) real, must be either purely real or purely 
imaginary. If « is real we can take it as positive without loss of generality, 
and equation (8.2) then shows that 6, must be zero, since Z must not tend 
to infinity as z> —«. But, if 5, is zero, (8.1) requires that a Au, 
which is impossible since Ay is essentially positive. Hence « cannot be real. 
We therefore take « = iw, where w is real; equation (8.2) can then be 
written in the form 
Z (C cos wz dD, sin wz)e Br (8.3) 
where C, and D, are real. The boundary condition (8.1) now requires 
D,w Au. (8.4) 


Hence the current density at any point inside the conductor is given by 


: Se inita Au . oP oP 
1 KE KC ¢ (cos we f sin w2}( > : 0). (8.5) 
Ww \CYy Cx 
where B c7(\?+-w") darcy. (8.6) 


Outside the conductor, equations (6.12) and (6.14) show that the electric 
field is . ‘ 
E CO. Bt ae CP 4 x (8.7) 
cy Ox 


The corresponding magnetic field is found from (6.15) and (6.16); these 
equations give on integrating with respect to ¢, 


cl, 3 ; P oP\ P 
nH = bes ain ewe t Au cos ew2)(* ; : ). (Au SIN wz wcosws)P' 
B | cx cy} w } 
for 2 6, (88) 
cC,A 2 ’ 
and H ——ePorad(e-“P) for z> 0. (8.9) 
p 


Equations (8.5) to (8.9) give the free modes of decay of the first type. 
In these modes the currents all flow parallel to the bounding plane. There 


are modes corresponding to all real values of A and w from 0 to ~, and to 
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il functions P(x, y) which satisfy (6.11), the distribution of current density 
n each mgde being a simple periodic function of z. 

The free modes of the second type are given by equation (7.8) and the 
elated equations of section 7. To obtain expressions for the current 


listribution and field in a form directly comparable with those given above 


for the modes of the first type, we write Cy B,w/A in (7.8), so that the 


rent density at any point in the conductor is given by 


;, ' 2 cP oP »P 
1 KE KC, ¢ “| COS w2 » COS Wz ; 
7 CX cy Ww 


sin wz}. (8.10) 


nd the electric field outside the conductor is given by 
E = (,egrad(e-“P). (8.11) 
The corresponding magnetic field inside the conductor is easily found 


be 
tarK - cP cP 
H C,e-P sin w2 ' 


WwW CY CL 


, 83. (8:59) 


hile outside the conductor there is no magnetic field. 

In these second type modes there is, in general, a component of the 
rent density in the direction normal to the boundary plane z = 0, 
ough this component is always zero at the boundary itself. There are, 
ist as in the first type, modes corresponding to all values of A and w 
between 0 and «, and to all functions P(x, y) satisfying (6.11). Also the 
urrent density is again a simple periodic function of z. 

There is, however, one peculiarity about these modes which requires 
explanation. The electric field outside the conductor, given by (8.11), 
rresponds to a distribution of charge on the surface z = 0, since by 
3.6) there can be no volume distribution of charge in the conductor. 


his surface charge distribution is given by 
o = E,/4n Cie PtP(x, y). (8.13) 


nd, since it varies with the time, it implies a normal component of current 
t the surface of the conductor of amount 
t Ms 
Cle P'P(x, y). (8.14) 
$77 — 
in apparent contradiction with (8.10), which makes 7, zero at z = 0. The 
alue of i, given by (8.14) is, however, of the same order of magnitude as 
the displacement current D/4z, which we have already assumed can be 
ieglected in comparison with i in the fundamental equation (3.1). Hence 
there is no real contradiction in the above results. Equation (8.10) is, in 


lact, an approximation, which holds good for the range of values of 8 (i.e. 
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of the conductivity «, linear dimensions A~!, w~!, ete.) contemplated, but 
which will break down if the value of £ is too large. 

The physical interpretation is that a minute portion (only) of the (posi 
tive or negative) current, flowing towards the boundary in any one of these 
modes, will carry on right up to the boundary, and there build up a surface 
charge, whose electric field is sufficient to deflect the remaining, and much 
greater, part of the current away from the boundary, so that at the 
boundary itself the current flow is practically parallel to it. The electric 
field E due to the surface charge o will be 470/e, whereas the magnetic 
field H of the corresponding current ¢ is given by curlH trac. Hence 
although the electric field of the surface charge is of fundamental impor 
tance in determining the distribution of current in the conductor. the 
magnetic field of the current associated with this varying charge is 
negligible. 

It is of interest to note that the free modes corresponding to the two 
different types of solution become identical when A = 0. Thus, taking the 
first type of solution and putting A = 0 in (8.5), we obtain 
By oP, oP, 


1 KC’, e- Pl cos wz 


, 04, (8.15) 


cy Cr 
2p 2p 
OP, oP, 


ox" cy* 


where (0. (8.16) 


so that P, is the real part of a function, P,+iP, say, of «+iy. Hence 


oP, cy oP, cx and éP, ox eP,/cy, so that (8.15) can be written 
: : oF. oP, 
1 K( ié BY cos eo2| =. = OO}, 
Ox cy 


and this is obviously of the form assumed by the second type of solution 
(8.10) when A = 0. In these particular solutions there is no magnetic field 
outside the conductor, and the electric field outside is everywhere parallel 
to the surface and independent of z. 

We may also add that the solutions (8.5) give the free modes of decay 
of current in a uniform plate of thickness h, provided w is now restricted 


to the discrete set of values nz h, where n is any integer. 


9. Periodic inducing field 

If the external magnetic system is periodic in the time, it may be 
analysed into simple harmonic constituents, and the effects of the separate 
constituents may be superposed in virtue of the linearity of the equations. 
Taking, therefore, a single harmonic, we now calculate that part of the 
induced field which has the same harmonic time factor. The remaining 


part of the induced field will be a transient part depending on how the 
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external variations were first set up. This part depends on the free modes 


fdecay of current in the conductor: it is negligible if the periodic external 


system has been in existence a sufficiently long time. 


If the period of the inducing system is 27/p, the equation (6.19) for the 
otential of the total magnetic field outside the conductor may be written 
nthe form 
ie Q Ae + Re-*)eit P(x, y), (9.1) 
here 4d and B are complex constants and the real part of the expression 


nthe right is to be taken. The modulus of the ratio 6/A will then give 


the ratio of the amplitude of the induced magnetic field to that of the 


nducing fields, and the argument of B/A will give the phase difference 
f the fields 


Since E must also contain the time-factor e’// the equation (6.6) reduces 


CZ / torKkptp 


~—| —=—- t-A2}Z = 0, (9.2) 
ga 
he solution of which is 
VA (ae® + he-")etvt, (9.3) 
er ( torkpip c*+°, (9.4) 
l ; a , ° - 
so that d 1$(x#+-A4 A?! (a4 | \4) A?! |, (9.5) 
v2 
here y~ firK Lp Cc". (9.6) 


In (9.3) we can obviously take the real part of @ to be positive without 
ss Of generality. The boundary condition (6.18) then shows that 6 must 
be zero, and (6.14) and (6.17) together with (6.20) require that 


cAa ip(A B). cOa pip(A | B), 
so that oo 
a ~ a Mo (9.7) 
c(@- Au) 
J 4 
ba oe. (9.8) 
0 Au 


Thus the scalar potential of the induced magnetic field outside the con 
luetor is given by 
6—Ap 


0 Au 


Le -A=¢'! Pla, y), (9.9) 


vhile the electric field inside the conductor is given by 


) 2 , 
E “HAP Acivt 06, of 0), (9.10) 


cy Ox 


c(é Au) 
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The current density at any point inside the conductor is given by 


. . y 2 \p 
i «E 2Kpip Aetvt a( m... ol. 


(9.11 
c(6+Ap) cy ox 


10. Periodic linear current parallel to surface of conductor 
Suppose that the inducing field is produced by a periodic current 
cece” e.s.u. flowing in the direction of Ox along the line y = 0, z — h. 
The magnetic potential Q, of this current is 2c, e'” tan-{(h—z)/y}, and 
it is easily found by Fourier’s integral theorem that this function may be 
expressed in the form 


OQ, 2 | cyexp(ipt+aAz Ah)sin dy, when z<h. (10.1) 


0 
Comparing this expression with (9.1), we see that P(a,y) is now sin Ay 
and A is 2cge- dd/A so that, from (9.9), the potential Q, of the induced 
tield is 


I~ 
be 


°@ D 
, ’ | 64 WHegexplipt Az—Ah)sin Ay - (10.2) 


Also, from (9.11), the induced current-density is given by 


. oe 
. 9 | aKpUp 


*») exp(tpt—Ah-+ 0z)eos Ay da, 10.: 
“(04 ru) xp(27 2 )cos Ay ¢ (10.3) 


0 
with 7, and 7, both zero. 
10.1. The induced field 


The components of the induced field are 


a, =9 
AQ °@ : 
H, : WC | Ae ex p(ipt—Az—Ah)cos Ay da, 
cy | O+Ap , 
QQ. , 6 4 ° . 
H. “ 7 ae | e Weexplipt Az—Ah)sin Ay dd. 
0 


In these expressions @ is a complex number and the real part of the 
integrand is to be taken. We can write* 


a {0 Au ° . 
H,--iH 2} re f cyexp(ipt) exp! Nh+z+iy)' da. (10.4) 
; : \0-+-Au 

0 
* Note that the direct field of the linear current can similarly be expressed in the form 
H iH, 2¢) COs pl (h ra “y). 
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Taking c, real (so that the inducing current is c, cos pt), and assuming 


the magnetic permeability , to be unity, we find that 


§@—Au | 27d? 21r6 . 
sl c,exp(ipt)' cyre (1 — | - exp(ipt)! 
\O-+Aw | 2 y? | 
C, COs pet | uv24(1 ut) w=} | | ey sin pt{ 2a? uv2i(1 T ut)! T uP) 


here u A/x, and « is given by (9.6). 


Hence the expression (10.4) becomes 


0 


H +-iH 2c, cos pt | |l—uv2{(1+-ut)!—au?h les" du 


2cyasin pt | [2u2—wv2{(1-+-u!)!+-u2} Jeo" du, (10.5) 
vhere C v(h-+-z)+-iay = &+1n, say. (10.6) 
Asymptotic expansions of the above integrals may be found by expand- 
ng the coefficient of the exponential factor in each integrand in ascending 
powers of uw, and integrating term by term. This leads to 
H +iH, = 2cya cos pt{l-1— V2(¢-2—3£-44 15-4 8156-8 + ...)} 


14.15f-6_315f-8+...)}. (10.7) 


2c, «sin pt{4f-3— 1 2(f-? + 3 


ae 

It will be observed that, when « (and therefore «) tends to infinity, 
H,+iH, tends to (2c, cos pt)/(h+z+ ty), so that the induced field is the 
same as that of a current, equal and opposite to the inducing current, and 
flowing along the line y = 0, z h. This is in agreement with Maxweil’s 
image-theory* for a plane conducting sheet. 

The expansions obtained above are only of value when |¢| is greater 
than about 6 if results to three significant figures are required. When ¢ 
is less than this it appears necessary to evaluate the integrals numerically. 
For moderate values of «, it is precisely these smaller values of ¢ for which 


the calculations need to be made, since the induced field is then only 


important for sufficiently small values of h+-z and y. 
If we write , ae : 
H, 2Cy {24 (E, cos pt t Y3(€, n)sin pt}, (10.8) 
H 2¢) a{Z,(€, n)cos pt+ Z,(€, n)sin pt}, (10.9) 


ve find easily from (10.5) that 


Y5(€, 7) | [l—wv2f(1+-ut)t— ut |e "cos un du, (10.10) 


Y(E. 7 [ 2u2?— wr 25(1+- ut)? + 2! le "cos un du, (10.11) 
ind Z,(€, »), Z.(€, 7) are similar integrals with cos wy replaced by sin wy. 


Clerk Maxwell, Proc. Roy. Soc. 20 (1872), 160-8. 
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[t will be observed that a know ledge of the above four 
two independent variables (€, 


functions of the 
7) is sufficient to determine the field at any 
point (y,z) for any set of values of x. p. and h. 

The four functions Y,. ¥2, Z,, and Z, were calculated for a series of values 
of (€,) from (0,0) to (10,10) and are shown in Table [. 


TABLE | 


Values of Y,. Y,, Z,. and Z, 








n 
I ; . 5 I 
: art } 105 I O1Q5 11 O5 14 
0 -_) O°47 o*1090 I O°O7 5 | 12 
Z) ) 164 271 0°25 2 I 
Z> 197 235 18y 1tt 006 } 
1 "401 5 345 ¢ O1g4 O°125 ol 
=e Yo 25 I O'l4y }f 0°026 0°035 3 co} | 
s) 7 - >> i 0°20 - 
L, O°123 189 O°221 0*209 155 0) 
Zs I 145 “1 3¢ O'103 006 5 
1 324 7 305 19 orgs ohn meta) )4 24 
) 17 128 5 O° o0158 123 
1:0 7 
Ly 048 0°02 [43 O°175 1O¢ , 
; 
2 53 Q 3 Itt o'066 14 0 
, 5 49 245 I rz Ol 37 104 
2-0 Ys O°0Qg5 3 OdQ9 023 O’oo! It Ol 
Z1 0-058 orl O°1 26 O°1 32 I ( 
Ze } 46 57 “063 0'055 } 
yy 2 C 196 18 O'15 o"129 I fe) 
} 64 ( 0°059 Oo Oo! 0°008 3 7 
30 7 ‘ oo 
“4 "OW ‘O75 0°0Q2 oO"102 10¢ } 
7 r I 0 _ > ny 
Zs 0°0 2¢ I 038 O'o4! j 
1 17 168 165 157 140 O'121 103 
1 
4:0 “a 45 17 0°0.42 029 o°018 Oo10 } 
Ly 15 0028 52 067 0-080 0°087 
Z ¢ O12 24 "027 0°03 D0 2¢ 
3 146 145 140 [41 129 O'l12 098 0°0467 
5-0 2 32 32 0 3¢ O°O17 O°o12 ( 04 
1 020 0°0 36 0°050 0°065 ” 0707 
2 ) _¥) 020 o"o2!I Cc 
) 0862 Sé 0°0857 08 39 731 049 
Vv. ‘O10f oO1Os 0°0103 0:0007 < 00 
2 7 2 
10'0 7 zs ae 
“y } 73 O143 I }2 
a g 0017 0034 re oof 





For values of |é+ in| less than 6, the infinite integrals for these functions 
were evaluated by numerical integration. The actual range of integration 
taken was from u — 0 to u 10, and the ordinates were calculated to five 
significant figures for intervals of 0-1 from xu 0 to u 


u l 5, and of 1-0 from x 


1, of 0-25 from 


to u “i 5 to 


of 0-5 from u 


2 to u 
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10. The integrals from u = 0 to wu = 10 were then calculated by 
Simpson’s rule and rounded off to three figures. An estimate of the 
remainder of the integral, viz. the integral from 10 to «0, was obtained by 
noting that the expressions in square brackets under the integral signs in 
10.10) and (10.11) tend to the values 1/8u* and —1/4u?, respectively, 
vyhen « becomes large. For example, the remainder of the integral for 
/,(0,0) is approximately | (1/8u*) du = 0-00004; this is neglected in com- 

10 


parison with the integral from 0 to 10, the value of which is 0-471. 

For values of |€+-iy| greater than or equal to 5, the functions were 
alculated by means of the asymptotic series (10.7). This series gives 
results correct to three places of decimals (i.e. within about 1 per cent. of 
the actual value) when |€+-i7 5, and gives results correct to four places 
vhen |€+77 10. 

In those cases which were calculated by both methods (i.e. where |€-+-7y 
is about 5), the results agree to within 1 or 2 per cent. It is believed that 
the results given in Table I are accurate to the three figures given except 
when |€+-in 2, in which case the last figure may be doubtful to the 
extent of 2 or 3 

Values of the functions when € or 7 is greater than 10 can be very quickly 
obtained from the asymptotic series, which is then rapidly convergent. 

The graphs of Y,, ¥;, Z,, Z, as functions of », for different values of & 
are shown in Figs. 1 and 2. Since & (h+-z),/(4aKp) and » = y,/(4zKp), 
each curve may be regarded as showing the distribution of a component 
of the induced field over a certain plane parallel to the surface of the con- 
ductor, the scale depending on the values of c,h, and ,/(47«p). For example, 


if /(4np) = 2/h, the curves ¢ 2 give the distribution of the field over 
the surface z 0, the scale for y being y inh, while the curves & 3 
give the distribution over the plane z = 4h, with the same scale for y. 


In any given case, the values of h and ,/(47x«p) will fix the € which corre- 
sponds to the surface z — 0 and the scale for y. Each larger value of é will 
then give the distribution of the field over a plane parallel to and above 
the surface. The curves show how, on moving away from the conductor, 
the induced field becomes less intense in the immediate neighbourhood 
of the line current but is spread over a much wider area. 

10.2. T he induced current distribution 


The expression (10.3) for the induced current simplifies to 


(9—A)cy exp(ipt+0z—Ah)cos Ay dd, (10.12) 


when b l. 
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Fic. 1. Graphs showing the distribution of the tangential component of the 
induced field over planes parallel to z = 0. The component is given by 
iM, 4c, ./(aKp){Y,(€, n)cos pt ¥,(€, n)sin pt}, 


where & = 2(h+2),/(7Kp), n = 2y,/ (Kp). 
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Therefore L(é.») = | ¥(€.) dé, (10.14) 








INDUCTION 





IN A SEMI-INFINITE CONDUCTOR 












































Se Sm 


Values of Z, 





Fig. 2. Graphs showing the distribution of the normal component of the 


induced field over planes parallel to z 0. The component is given by 
H. $09 ./ (KP) Z,(€, n)cos pt Z,(E,n)sin pt}, 
where é 2(h +-2),/(arKP), n 2y,/(7KP). 
\t the surface (z 0) of the conductor, this becomes 
Cy a* 4 ‘an 
°— cos pt [{(1+-w*)?+- u?}t— uv2]e-%" cos nu du+- 
CTNZ 
0 
c y? ss r , 
= -sin pt | {(1- u*)t—yrite-%" cos nu du, (10.13) 
CTINZ 
0 


n substituting the value of @ and taking the real part. 
The integrals on the right can be computed in terms of the integrals 


i(é,m) and Y,(€, 7») already obtained. For, consider 


L(€, n) | [| 2u—~+2{(1+-w*)!+-u?} Je E" cos nu du. 
. 
Ve have 
el, eee ° , 2131,-£ , 
. [ 2u?—wv2{(1+-u*)!+ w?}4le-§" cos nu du Y3(€, 7). 
cet 
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since [,(€,n) > 0 as €-> a. Again consider 


I,(€, 7) ( V2{(1+-u*)!—u?} 


¢ 
*e-$" Cos nu du. 


Here 
ol, j - ms ¢ ree € 
c | uv 2}(1-+-ut)'—u*e-$" cos nu du Y,(€, 7) ego 
of , . rey 
0 
and therefore I,(€, 7) Fe -—Y,(é, ”) dé. (10.15) 
g SS TY 


This integral is readily proved convergent by writing Y,(€, 7) as the real 
part of the corresponding integral in (10.5). Incidentally we obtain 


Y,(€,) > a————_ i E+in| > 00. (10.16) 
ry 


The expression (10.13) for the induced current density at the surface 
may now be written in the form 


mi 


To determine the current density below the surface it is necessary to 
compute the integral on the right of (10.12). Apart from a factor «?, this 
integral is a function of the three variables ah, ay, and az. 


11. Some special inducing fields; indeterminate problems 
As noted in section 6, the fields which correspond to the case A = 0) 
require special consideration. We have already considered, at the end of 
section 8, the modes of free decay of current systems corresponding to this 
limiting case. We now consider the nature of the solution when there is 
a periodic magnetic field of corresponding form arising from external 
sources. The potential of such a field will be expressible in terms ot 
functions of the form 
Q, (Ayz+ B,)e'"'*P(x, y), (11.1) 
@P @P 


a2 !' 5,2 
ox" cy 


where 0. (11.2) 


This includes, as a special case, that of a uniform field normal to the 


surface z = 0; its potential is then —A,ze’”", Again, if the field is due to 
a uniform sheet current flowing above and parallel to the surface z = 0 in 


the x-direction, its potential is of the form — B, ye. 


i lf vh, xy)cos pt+-I,(ah, «y)sin pt}. (10.17) | 
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On putting A = 0 in (6.6) and assuming Z contains the periodic time 
factor e’?!, we obtain 
e?Z 4akpip 


> 4, 
CozZ* oC 
so that Z C exp{ipt+ (1+2)yz}, (11.3) 
where y \(27Kpp/c?), (11.4) 
the positive value of the square root being taken for y, because Z must 
not become infinite as z + —oo. Hence, from (6.13) the value of E inside 
the conductor is given by 
: ' oP oP - 
E Cexpizpt+(1- iva(T ——, OV}. (11.5) 
oy ox 
Outside the conductor E must be of the form 
. wer oP 
E (az Dyer ; in ~O}-7 (11.6) 
oy OX 


and the continuity of the tangential components of E at the boundary 
requires that C=b. (11.7) 
Using (6.15) the corresponding magnetic field inside the conductor is 

found to be 
cO(] OY exatins me ma(S. oP 0). (11.8) 


Ox Oy 


H 
Pp 
Outside the conductor the total magnetic field must have a potential of 


the form 


Q) (Az+ B)e'”'P(x, y), (11.9) 
so that H — eit] (Az py ot , (Az BP api. (11.10) 
Ox oy J 


The continuity of the normal component of B at z = 0 then shows that 
A = @, (11.11) 


and the continuity of the tangential components of H requires 


B ni—) (11.12) 
PP 

Now if we write the potential of the induced field as 
Q, (A,z+B,)e”"P(x, y), (11.13) 
we have A A,+A, and B B,+B,. Hence the above equations show 

thé 
” A, Ay, (11.14) 
and B, 1 — 7m. (11.15) 
PP 


From equation (11.14) we infer that, if the external inducing field given 
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by (11.1) contains a component normal to the surface of the conductor, 
this component will always be completely annulled by that of the induced 
field, which is equal and opposite to it. 

On the other hand, equation (11.15) is insufficient to determine the 
component of the induced field parallel to the surface, because it contains 
two unknown coefficients B, and C. Since no other necessary relation 
exists between these coefficients, we conclude that, apart from the above 
result relating to the normal component, the problem of finding the 
induced field and current system corresponding to an inducing field of 
the form (11.1) is not a determinate one. 

This result is obviously associated with the fact that the conductor 
extends to infinity, and the field at any point P is determined not only 
by the currents in the neighbourhood of P, but also by the currents in 
the outlying portions of the conductor. 


12. Conductor treated as a limiting case of a large spherical 
conductor 
Since any actual application of the solution of the present problem must 
necessarily relate to a restricted portion of the surface of some physically 
finite conductor, it is of interest to treat the conductor as the limit of a 
finite conductor of particular shape, when its size is increased indefinitely. 
It is sufficient to consider the case of a spherical conductor of radius a, 
which will subsequently be made to tend to infinity. In this case the 
potential of the external inducing field is conveniently expressed in terms 
of spherical harmonies of the form 
Q, -AetPlyng-"-lg , (12.1) 
where S,, is a surface harmonic of degree n. 
The potential of the corresponding part of the induced field is found to 
be (see, for example, chap. xxii of reference 6) 
Q, = — Beimly-n-rgn +38, (12.2) 


the coefficient B being given by 


2 ant 
os a | mua Bn ; (12.3) 
A n+] m(u—1)+-OF,_3(f)/F,44(C). 
where C2? = —4rxpipa?. (12.4) 
It follows from (12.3) that 
a en (12.5) 
A n+] 


Now the normal component of the inducing field at the surface is 
nAe'S, and that of the induced field is —(n+-1)Be'S,. Hence, in the 
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limit as @ 00 these components cancel each other. This is in agreement 
with the result of the last section. 

On the other hand, the tangential components of the induced and 
inducing fields are, ultimately, in the ratio n/(n+1). Now any surface 
harmonic of finite degree n will lead to a field which is sensibly uniform 
over a limited region of a sufficiently large spherical surface. If we then 
confine our attention to this limited region, we might be tempted to treat 
the conductor as a semi-infinite one, and the field as uniform. But the 
results of the last section show that, if we do this, the tangential com- 
ponents of the induced field become indeterminate. This is, of course, in 
agreement with the fact that the above limiting ratio for the tangential 
components of the induced and inducing fields in the case of the sphere 
is n/(n+-1), and n may be any positive integer, so that this ratio may have 
any value between } and 1. Moreover, by taking conductors of other 
shapes (e.g. a cylinder), it is easily seen that still other values for this ratio, 
outside the above range, can also be found. 

If in the above we remove the restriction that » remains finite, and 
allow » and a to tend to infinity together, we can obtain a field which 
varies sensibly over a region of the spherical surface small enough to be 
treated as plane Thus suppose 


nia>A as a->o. (12.6) 


Then if S, is the tesseral harmonic P%(cos @)cos(s¢é+e), we find, on con- 


sidering the limiting values of the expressions (12.1) and (12.2 


Q,/a > —Ae'eJ (A R)cos(sd +e) (12.7) 


and Q/a > Beie—*=J (XR)cos(sd+-e), (12.8) 


where the z-axis is taken along the line 6 = 0, and R is the distance from 
this axis. 

To obtain the limiting value of the ratio B/A given by (12.3), we make 
ise of the recurrence formula for Bessel functions, namely 


2n+-1 


J, _4(C) +S, .3(C) J, .4(G). (12.9) 


S 


Writing nB and denoting the ratio J, ,,(¢)/J,,,(¢) by p,, this formula 


nH” 


ves 


S I 2n+-1 
Pn1 , (12.10) 
p np 
Hence, if p, > p as n > &, we have 
| » 
p 5 
p B 


(12.11) 
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PRICE 


and the negative sign must be taken since p does not become infinite as 
8 +0. Hence from (12.3) we obtain 
B a Tn 2u 
A p—1+p%{1—./(1—B?)}" 
where @ is given by (9.4). 
This result is thus identical with the result (9.8), and affords a check 
on the previous work. 


6 Au 


] 
] OAAw’ 


_ y (1—f*)—p é 
JB) 


13. Aperiodic inducing fields 

When the inducing field has a time factor A(t), which is not simply 
periodic, the induced field and currents can be determined only if the 
initial conditions of the field are given. Since, however, the differential 
equations are linear with respect to ¢, the effects of different time variations 
of the inducing field can be superposed, each such variation merely adding 
a new current distribution to that already existing (though continually 
decaying) in the conductor. It is therefore sufficient to take, as the initial 
condition, that the current density is everywhere zero in the conductor. 
The induced field and current distribution, corresponding to any given A(t), 
are then determined by finding the particular combination of free modes 
of the first type (section 6) which, when combined with the particular 
integral of the equations corresponding to A(é), satisfies this initial con- 
dition. This can be done with the aid of the Fourier—Bessel integral, but 
a more convenient method is to use the operational calculus. 


14. The operational method of solution 

A formal solution of the differential equations and boundary conditions 
is first obtained, treating the operator D = 0/ét as an algebraic constant. 
Thus solving equation (6.6) in this way, we obtain 


Z(z,t) = e%a(t)-+-e-"b(t), (14.1) 
where 62 = 4rxpD-+d?. (14.2) 
To satisfy the boundary conditions at z = 0 and z = —oo, we find, 


exactly as in section 9, that 5(¢) must be zero and 


2uD ‘ 
a(t . A(t), 14.3) 
at) = — Al 
0—Au 
Bit A(t). 14.4 
(t) 0 Ay (¢) ( ) 


The coefficients of A(t) in the above two equations may be regarded as 
operators, which, acting upon A(é), give the required functions a(t) and 
Bit). The result of any such operation, say 


a(t) = F(D)A(8), (14.5) 
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nite as ] 4 
is given by a(t) 7 A(t—r) f(r) dr, (14.6) 
0 
where f(#) is the solution of the integral equation 
F(D) = D | f(te-”* dt. (14.7) 
check Fs 
Thus a(t) and B(t) can be determined from (14.3) and (14.4) without further 
reference to the differential equations or boundary conditions. The value 
imply of Z(z,t) can then be obtained from (14.1) as the result of the operator e* 
if th acting on a(t). The induced field and current distribution are thus obtained 
- 
ential explicitl} as functions of the time. . . 
om The function f(t), which is defined by (14.7), is the value of a(t) for the 
ddine ease when the time factor A(t) of the inducing field is Heaviside’s unit 
seal function H(t), i.e. for the case when the potential of the inducing field 
nitial increases suddenly by one unit at time ¢ = 0. The equation (14.6) is the 
sien formal expression of the superposition principle. 
| A(t) 15. Evaluation of the induced field 
nodes From (14.4) the operator F(D), which is required to obtain the time 
cular factor B(t) of the induced field, is (@—Ay)/(@+-Au), where @ is given by 
con 14.2). Taking the case when pu 1, and expanding F(D) by the binomial 
I, but theorem in powers of D-', we have 
4 2A2 2AD (—)*-1(2n—3)!/ A \* 
F(D) l 4-3 , om | 
t7kD (47D) | 2"°n!(n—2)! \4anD | 
tions ; (15.1) 
ah Hence, using the known result 
tant . ' 
D-" H(t) t”/T'(m-+-1), (15.2) 
we find 
14.1) 2r7t 2, t (—1)"-1 A*t \" : 
i = 14 i }1 += ATV ..4. (15.3) 
14.2) dork «=o MN \K/ | n! (2n—1)(2n+-1)\47K J 
find. This series is convergent for all values of t. We can also write the result 
in the form g 
F(t) | } wD) (wu) (15.4) 
14.3) +} 
where D (uw) is Weber’s parabolic cylinder function, and 
{ Mt - 
14.4) u / . (15.5) 
A) \27K 
d as Using the asymptotic expansion for D,(u) we can express f(t) in the 
and 


alternative form 






; “ . 2 | 4! (2n+-2)! , — 
. f(t)~ u~% I —— +... (—)” ——+t 3, (15.6) 
14.5) NY (-} : : J 
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The two expressions (15.3) and (15.6) serve to evaluate f(t) conveniently 
for all values of f. 


16. Aperiodic linear current parallel to the surface 


In the case of a line current of strength C(t) along the line y = 0, z = h. 
we can express the potential of its magnetic field in the region z < h, as 


in section 10, in the form 


Q, 2 ( C'(t)exp(Az—Ah)sin Ay dX/A. (16.1) 


Hence the potential of the corresponding induced field is 


= 
: ME Cinexp\ Az—Ah)sin Ay dd/A, (16.2) 


I~ 
~ 
bo 


J O- Au 


and the components of the induced field may be written (ef. section 10) 
in the form 


H+iH, = 2 


u 


* 9) 
MH Cex! A(h+z+1iy)! da. (16.3) 
0 + Nu 


0 
To determine the effect of a sudden increase of amount Cy in the line 


current, we write C(t) C, Hit). (16.4) 


and we then have from (15.4), taking pu l, 


H +1 SC (4 ; D_.(u)ex | (Pra |2+ay)eu 1 2\ du. (16.5) 
/ 2 o,/ v | 3 Pia | d y) al 


If the power series (15.3) is inserted for exp(— u?/4)D_,(w) in this integral, 


we obtain an asymptotic expansion for 1/,+-7H_, valid for large values of 


(h+z+iy)t. This gives a convenient method of calculating the field for 
small values of ¢. For large values of ¢, it appears necessary to use the 
asymptotic expansion (15.6) to evaluate the integrand, and then to com- 
pute the integral numerically. 
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ON SOME DUAL INTEGRAL EQUATIONS OCCURRING 
IN POTENTIAL PROBLEMS WITH AXIAL SYMMETRY 


h 
hos | By C.J. TRANTER ((Wilitary College of Science, Shrivenham) 
Received 9 May 1950] 
16.1] 
SUMMARY 
The use of a Hankel transform can reduce the solution of Laplace’s equation in 
lindrical coordinates (p,z) in the region 0 < p 10, O Z h when the boundary 
16.9) ndition on 2 0 is a ‘mixed’ one and that on z h is of the usual type to the 
lution of the dual integral equations 
1 10) F ' 
| |G u)f(ued, pu) du q(p) (0 p l), 
0 
16.3) 
f(u)Jo( pu) lu i) (p 1), 
} 0 
line | where G(u), g(p) are given functions of the variables indicated and f(u) is to be found. 
16.4) \ formal solution of these equations is given and, as an example, the solution is 
ipplied to find the potential due to a circular disk at constant potential placed with 
s plane parallel to, and equidistant from, two earthed parallel plates. 
16.5) : 
1. Introduction 
PoTENTIAL problems with axial symmetry for an infinite plate involve the 
gral, : ae ‘ms : ; ' 
solution of Laplace's equation in cylindrical coordinates (p, z) in the region 
Ss ot ° . b> 
: 0<p 0,0 <2 < h (the plate thickness) under given conditions on the 
for eas 
: plane boundaries z = 0, z = h. If the boundary conditions are such that 
the 


either the potential, its normal derivative, or a linear combination of the 


‘om- 


potential and its derivative are specified over the whole of each plane face 

of the plate, a convenient method of solution is by the use of the Hankel 

transform. In such cases the use of this transform eliminates the radial 

variable and the problem reduces to the solution of a simple differential 

equation in one independent variable z with known conditions on z = 0, 
h. 

If, however, the condition on one boundary is a mixed one, for example, 
if the potential on z = 0 is a prescribed function of p inside a certain circle 
and its normal derivative is zero outside this circle, the solution by a 
Hankel transform reduces to that of a pair of dual integral equations. 
\ well-known example of this type is the field due to a single circular disk 
at unit potential. Titchmarsh (1) uses the Hankel transform to reduce this 


(Quart. Journ. Mech. and Applied Math., Vol. 3, Pt. 4 (1950)] 
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' hs ' 
problem to the solution of the dual integral equationst 
of | 
uf(u)Jy(pu) du lL (Q0<p<)l), of 
9 (1) | 
flu )J( pw) du 0 lo > -B). 
0 
In this example, if V is the potential, the boundary conditions are 
V lL, @- <a, 
p when z 0, (2) W 
OV /ez 0, p l 
LO 
and V>0O as z2>0. (3) q 
For the similar problem of the field due to a disk with its plane parallel to, | t 
and equidistant from, two parallel plates at zero potential, the boundary | & 
conditions are (2) and, instead of (3), | ir 
, V 
} 0 when z=h, (4, ° 
re 
2h being the distance between the plates. This problem is considered in | _,, 
section 3 and is there reduced to the solution of the dual equations | 
c | 2 
| u-! tanh hu f(w)Jy( pu) du 1 (0<p<l), | 
0 (5) | 
| f(w)Jo(pu) du 0 (p> 1). 
0 
The most general integral equations of this type for which a solution \ 
is available are, as far as I am aware, the pair considered by Busbridge (2) |‘ 
‘ , i 
| uf(u)J,(pu)du=g(p) (0< p< 1), 
0 (6) | 
| f(u)J, (pu) du 0 (» > I}, 


0 
in which g(p) is a given function of p and a, v are suitably restricted. These 
equations are not sufficiently general to include the pair (5) which is typical 
of potential problems in a region of finite thickness with one ‘mixed’ 
boundary condition. 


+ The equations actually derived by Titchmarsh are 
s 
F(u)J,( pu) du 1 (0 p ys 


L 


u F(u)J( pu) du 0 (p> 1), 


which can clearly be put into the form (1) by writing f(w) 


u F(u). 
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For potential problems with axial symmetry, only the Bessel function 
of zero order will be involved, and I have attempted in section 2 a solution 


f the dual equations 


| G(u)f(u)Jo(pu) du = g(p) (9 <p <)), (7) 
0 
| f(u)J,(pu) du 0 (p £, (8) 
0 
where G(u), g(p) are given functions of the variables indicated and f(w) is 
to be found. The solution given is a formal one only and the difficult 
juestion of convergence is not considered. Under certain circumstances 
this solution is likely to be suitable for computation purposes. As an 
example, the potential problem involving the solution of Laplace’s equation 


n cylindrical coordinates under the boundary conditions (2) and (4) is 


worked out in fair detail in section 3. That the method of solution of 
equations (7) and (8) could be extended to cover the case when J, (pu) is 


replaced by J,(pw) seems fairly ciear but is not attempted here. 


2. The solution of the dual integral equations (7) and (8) 
We start by considering the integral 


I(m, k, p) - KJ, (U)Io( pu) du, (9) 


m 


where m is a positive integer (or zero) and k is real and positive. Under 
these conditions the integral converges for both 0 < p l and p > land 
ts value is given by Watson (3a). When p > 1, there is a term ['(—m) in 
the denominator of the factor multiplying the hypergeometric function in 
Watson’s formula and hence the integral vanishes for such values of p. 
By taking 

f(u) yi-t > An JFom+ihU)s (10) 


ind by the use of the result just stated, we see that the second of the dual 
equations (8) is satisfied by this choice of f(u). 

The coefficients a, have now to be chosen so that the form of f(x) 
wssumed in (10) satisfies the first equation (7). For this purpose we need 
a preliminary result and we first consider the value of J(m,k,p) for 


0<p 1. Again quoting from Watson (3a), we have 


m+) EF (m+1,1—k—m;1;p?) (0 < p <1). 


I(m,k, p) Ik-1P'(m ti 


The hypergeometric function can be expressed in terms of a Jacobi 
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polynomial by first using the transformation (see, for instance, Magnus 
and Oberhettinger, 4a) 
F(a, b;¢;z) = (1—z)*-*~", F,(c—a, c—b; €; z). 
This gives 
T(m+1) 


I(m, k, p) 2k-AP(m 4 b" p?)k¥1,F(—m,m+k:1;p?) (0 <p <1). 


Introducing Jacobi’s polynomial defined by 


F, (x,y, x) = .Fj(—m,a+m;y;2) 
(Magnus and Oberhettinger, 45), this can be written 
C(m-+-1 i 
I(m, k, p) ei ) (1 p?)/ FZ (k, 1, p”) (0 <= p <= 1).. 


2k-1P'(m-+k) 


By applying Hankel’s inversion formula to (9), using (11) and the fact that 


I(m,k,p) = 0 for p > 1, we have, as the preliminary result required, 
r 
. (m T 1) % O\L az 9 6 
u Jom iK(U) Dk 1D (m + k) p(l p*)' LA Ak, 1, p*)Jq(up) dp. (12) 
0 


Substituting from (10) in (7), we see that a,, has to be chosen to satisfy 


x 


> e ( G(u)jul-*J,,,, .;.(u)Jg(pu) du = g(p) (0 <p <1). 
m 0 * ‘ 
0 
Multiplying by p(1—p*)*-".4,(k, 1, p?), where v is a positive integer or zero, 
and integrating with respect to p between 0 and | and assuming that the 
order of integration can be interchanged we have 
. 1 


= ¢. ( G(u)u* Jy, .,.(u) du ( p(1—p?)! LF (kh, 1, p®)Jy(up) dp 
m=0 5 i 
: 
| 9(p)p(1—p?)*-1.F,(k, 1, p?) dp. (13) 
Writing, for brevity, : 
: 1 
— (n+1) fF nb .8 sk P 
E(n, k) KAP(n w | g(p)p(1— p?)*-LF,(k, 1, p?) dp, (14) 
0 
and using (12), (13) gives 
> Qn [ G(u)u' lis i(U)As,, ile) du E(n, k), (15) 
m=0 : 


0 
which is a set of simultaneous equations to determine the coefficients a,,. 
Since (Watson, 35) 


u Jom (Ug, n(u) du 


0 
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equations (15) give simple expressions in closed form for a,, when 


9 


G(u) yrk—2_ 
part from being real and positive, k is at our disposal and the solution 
f equations (15) is facilitated by trying to choose k so that the difference 
(u)—u?*-? is fairly small. For the type of application we have in mind 
this is generally possible, for G(uw) will usually contain, besides u, a para 
neter h and will behave like a power of u for a certain value of h. For 
example, in the pair of equations (5), G(u) = u-! tanh hu and G(u) behaves 
ke u-! for large h. If then we are interested in the solution of (5) for 
wee h, we take k } and rearrange (15) so that the difference G(u)—u7! 
sinvolved. Returning to the general case, we write (15) in the form 


(U J, y\U) du 


} 


“= 2h G'(u) 1 tau =o (Ue, 


2m 2n 


,(u) du E(n,k). (17) 


Vriting 
Ln, (4n+-2hk) ( fu2-2*G(u)— Liao, 4 4.(U)Jo,, .,(u) du, (18) 
0 
nd using (16), we have 
An,+ > Land (4n+-2k)H(n, k). (19) 


By a proper choice of & the values of Z,,,, can be kept fairly small, and 


the coefficients a,, can be found from the set of equations (19) by an 


terative process. Expressions for the a,, calculated in this way are 
necessarily somewhat complicated in the general case and it does not 
seem worth while to give them here. That the method outlined can give 

useful solution is illustrated in the example considered in the next 


section. 


3, An example 

Here we find a solution for the potential due to a circular disk at unit 
potential placed with its plane parallel to and equidistant from two 
earthed infinite parallel plates. The form of the solution is well suited to 
lumerical computation when the parallel plates are fairly far apart. 


We use cylindrical coordinates (p,z), take the disk as z = 0, p < 1 and 
the plates as 2 Lh. If V is the potential we have to satisfy 
a] 1 eV ay 
; : — = @ (20) 
Cp" p cp Coz* 
vith 


when 
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since there is symmetry about the plane z = 0, and 
V=0, O0<p<o, when 2<=A. (22 
The same solution applies to the steady temperature in an infinite plate 
of thickness h when its lower surface is maintained at zero temperature 
while, on its upper surface, the temperature is kept at unity inside a circle 
of unit radius and there is no loss of heat outside this circle. 
Let V be the Hankel transform of V, i.e. 


V ( pV J(up) dp, (23) 
0 
then 
o2V > @yV -{ &V av 
; | p-— Jy(up) dp | (p s+: }Ja(up) dp. 
Cn" R cn" ‘ cp" Pp, 
0 0 
using (20); also 
i eV i Cc V , ) 
| Pa Jy(up) dp | (Jo(up) I upd (up); dp, 
e cp” J 2p 
0 0 
integrating by parts; hence 
aad | "ov , ae - ' 
~ | upd (up) dp U | J {J (up) | upd o(up)} dp, 
oz" J of P 
0 0 


integrating by parts again. Using Bessel’s equation, the expression in { } 


is 


in the last integral is —wpJ,(up) and we have 
eV /éz2 = wd. 
The solution of this equation is 
V Ae**+ Be-, 
and since from (22), V 0 when z h, A -Be-2"h, Hence 
V Bfe-¥=— eue—2n)\ 


The inversion formula for (23) gives 


V ( uBle-"=— eG-2 J ( ow) du, (24) 


0 
and substitution in the boundary condition on z = 0, equations (21), give 


[ uBi1—eh"\ J,( pu) du | @<p< Bj, 
0 
| we Bil+e*"|Jy(pu)du = 0 (p> 1). 


0 
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u> BSL 


e—2h os 


we have for solution the dual integral equations 


( ul tanh hu f(u)J,( pu) du 1 (O<p<l), 


0 


[ T(u \J,(pu) du 0 (p 1). 


0 


These are a special case of the more general equations already discussed 


1 which G(w) 


tends to infinity, and we are interested in a solution for large h, we take 


From (10), (19), (18), and (14) the solution is given by 


J(u) u > On Jomsi(X), 


m= 


vhere the coefficients a,, are to be calculated from 


A+ YX Linn dm = (4n+1)E(n, }), 


? 0 
which 
L = (4n+-1) | {tanhhu—l}e MTom sy(U)Joy .4(u) du, 
v0 
vd 1 
E(n,} 221 (n+1)/T'(n 5)} p(1 p”) !F (4,1, p*) dp. 
0 
Now (see Magnus and Oberhettinger, 4a), F(a, y, x) l and 


l 
ia r)°- VF (a,y,2)F(a,y,2%)dx=O0 (m ~n), 
0 
) that H(n, 4) vanishes for n 0, while for n 0 


1 
E(0, 4) (2/7)? | p(1—p?)-* dp (2/7r)?. 


m=0 
The iterative solution of this set is 


’ a [ ¥ LA va a 
a (</7r)*| 0 C,+C,—Cr+...], 


u-'tanhhu and g(p) = 1; as G(u) tends to u-! when h 


(25) 


(26) 


(28) 


(29) 


(30) 


(31) 


(32) 











418 C. J. TRANTER 


where 5,, = 1 for n = 0, 6, = 0 for n > 0 and 


7 


n Lon | 

( m > Ri ( m? 
— 4 (33) 

( n : 3 | ( m? 

m=0 | 

; \ 

J 


and so on. 


A convenient way of evaluating L,, ,, is to obtain it as a series of powers 

of 1/h. Writing the hyperbolic function in (29) as a series of exponentials, 

we have 

(8n+2) > (—1)* | e-""J,,, 3 (u)do 
s8=1 


I 


“mn ~ 


ny(uju-? du. 
0 


The infinite integral can be expressed as a power series in 1/(2sh) by using 


a result given by Watson (3c). If we do this, interchange the order of 


summation and use the result (Whittaker and Watson (5)) 


> (—1)-1s-” = (1-2! *)(p) (p > 1), 
s=1 
log, 2 (p = 1), 
where ¢(p) is Riemann’s zeta function, we find 
(8n+2) ¥ (—1)""1> 
r=0 
I'(2m+ 2n+-2r+ 1) (2m + 2n-+ 2r+ 2)(1 — 2-2m—2n-2r) 


» Dam+an+4r+2(r) I D(2m+-r+ 3)P(2n+7r-+ 3)P(2m+ 2n+r- 2) 


“m,n 


—C(2m+-2n-+ 2r+-1) 


h2m +2n+2r+1 


(34) 


For Ly the first term of the series involves the limit as p approaches unity 
of (1—2!-?)(p); this, of course, must be replaced by log, 2. 
If h is so large that powers of 1/h above the third may be neglected we 
find from (34) 
37 Log = (log, 2)/h—£(3)/(8h3), dirDo, = €(3)/(8h), 


— rly. = &(3)/(40h), 


and that all other L’s are negligible. Substitution in (33) then gives 


Co = (2 7)| (log, 2)/h—Z(3) (8h3)], 
C (2/)¢(3)/(8h?), 


Co = (4/2?)(log, 2)?/h?, 


3 /7)(log, 2)3/h3, 
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ON SOME DUAL INTEGRAL EQUATIONS 


and all other C’s negligible to this order. Equations (32) then give 
Ay (2 7)*{ I (2 log, 2)/(ah) 


(2 log, 2)?/(ah)?+-(2 log, 2)3/(a#h)’— 2C(3)/(8h)], 


To this order, (27) gives 
f(w) wag J(u) a, J;(u)}. 
B can now be found from (25) and the potential itself is then given by (24). 
Ifh is such that higher-order terms must be retained the algebra becomes 
somewhat complicated and it would probably be better to draw up tables 
of the Z’s and C’s as numerical functions of h. 
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NOTE ON DIFFRACTION BY AN EDGE 
By D. 8. JONES (Dept. of Mathematics, The University, Manchester) 
[Received 6 June 1950] 

SUMMARY 


On the assumption that the electric current and charge on a perfectly conducting 
sheet are integrable in the neighbourhood of the edge, it is shown that the line 
distributions of electric current and charge on the edge are determined by the 
boundary conditions on the sheet. As a boundary condition on the edge it is postu- 
lated that there be no line distributions of magnetic current and charge; it is shown 
that the field is then determined uniquely. It is also proved that, if the currents 
near the edge can be expanded in power series, the index of the first term is (2p+1)/2, 
where p is an integer, greater than —2 or —1 according as the current parallel or 
perpendicular to the edge is considered. Finally the diffraction, in three dimensions, 
of a plane wave by a semi-infinite plane is obtained. 


1. The expressions for the field 

Ir has been pointed out by Bouwkamp (1) and Rayleigh (2) that the 
problem of diffraction by an edge is not solved uniquely unless some 
boundary condition on the edge is specified so that the singularities there 
can be determined. In what follows the singularities formed by electric 
current and charge, magnetic current and charge are discussed. 

Let S be an unclosed surface bounded by a simple closed curve C which 
has neither double points nor discontinuous changes in direction. Let 8, 
be the surface generated by the circumference of a circle of small radius 6, 
whose centre moves along C. Let S, be that part of S which is not 
inside S,. We take S to be perfectly conducting and consider only periodic 
phenomena varying as e’, this factor being omitted throughout. The 
m.k.s. system of units is used. 

Let a field in which the electric intensity and magnetic flux density are 
E,, B, respectively fall on S, which is surrounded by free space. Let E, B 
be the electric intensity and magnetic flux density of the resulting field. 
Assuming that E—E,, B—B, satisfy the radiation conditions at infinity* 


tangential derivatives on S,+S, by means of the formulae 
— _ lw 
E = E,+ - 
4ork? 
| {—22yB, AdS+(dS. curl’ B, grad’ p+ (dS A E,) A grad’ d 
. | LB, . cur ,grac ——s 1) A gre , 
S ; S. 


* These conditions are not essential in sections 1—4 since E,, B, can be replaced by any 
ficld which remains finite at the edge without affecting the argument. 


[Quart. Journ. Mech. and Applied Math., Vol. 3, Pt. 4 (1950)] 





we can express E, B in terms of their tangential components and their 
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B= B,4 — " = Ey dS)ys—(dS \ B,) A grad’ (dS. B,)grad’y}, 


where k = 27/(wave-length), 4 = e-'*"/r, and 


r,y,2) and (a, y’,2’) being the point of observation and a point of S, (or S,) 
respectively (see Silver, 3). The primes on the vector operators indicate 
operations with respect to (x’,y’,z’). On S,, E, and B, are the electric 
intensity and magnetic flux density. On S,, nA E, and n.B, vanish, 
whereas in other terms B, is the discontinuity of the magnetic flux 
density across S,; dS (= ndS) is along the normal, out of the space 
under consideration, to S,; on S,, dS is along the normal, its direction 
being determined by the choice of B,. We suppose that the integral over 
‘, tends to a limit as 5 > 0, i.e. we take the electric current and charge 
on S, to be integrable in the neighbourhood of C. We shall in fact assume 
them to be Lebesgue-integrable. We further suppose that the integral over 
‘, tends to a finite limit as 6 > 0. This is certainly true when the fields are 
produced by line distributions of electric current and charge, magnetic 
current and charge on C’. In general there will be distributions on C of 
electric current and charge included in the integral over S, when 6 > 0. 
We consider the possible existence of line distributions (of finite magnitude) 


in excess of these. 


2. The line distribution of electric current 

If we consider the field of the line distribution at a point of S near C 
the dominant portion comes from the nearby part of C. From this point 
of view the most significant line distribution for the component of magnetic 
flux density normal to S is that of electric current since the nearby parts 
of the other line distributions effectively provide no contribution. Thus, 
supposing for the moment that only an extra line distribution of electric 
current is present, we have 

B B, (dS \ B,) A grad’ %-+-curl Ib ds, 
S C 

where / is proportional to the strength of the line distribution and the 
integral along C' is in the positive direction obtained by the right-hand 
screw law from dS. On S, but not necessarily on C, the normal component 
of B vanishes and this condition supplies a certain integral equation. 

Let O be a point on C and Ox that curve which is directed into S and 
which is the intersection of S and the plane normal to C at O. We consider 


the boundary condition at the point P, which is on Oz and is close to O 
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so that OP (= x) is smal]. The contributions to B at P from B, and those 
elements of S and C which are more than a certain small distance from 0 
will be virtually independent of x and will be o (1/x) as a +0. For the 
contribution of the near elements we shall obtain a sufficient approxima- 
tion if we take x so small that the near elements of C can be regarded as 
lying on a straight line of length 2m with O as mid-point, while the near 
elements of S lie in a plane rectangle of sides 27, «. Here 7, € are small, 
but are chosen so that 7 >« >a. The y-axis at O is taken parallel to n and 
the z-axis is chosen so that the axes form a right-handed orthogonal system. 

Let B = (B,, B,, B,), By, = (B,,, By,, B,,) and [,(2’,z’) = iwB,,/4ak?, 


; 


Lia’ 2°) iwB,,/47k* so that I,, 1, are proportional to the electric 


1x? 


currents in the near elements of S. The contribution to B,, at P of the 
near elements is 


€ 7) 7 
oO. Thess i see ee ee ee 
lim - | | IT (« ,2)——~—ZLi{2’,z’) | dz'dz’ —— T(z’) dz’, 
y—0 WW . | C2 Ox | Ox , 
0 n 7 
where y’ = 0,2 = Oin &. We assume that ol,/éx’, L, are continuous and 


differentiable except possibly at x’ = 0. We further assume that ¢J,/é2’, 1. 
and their derivatives with respect to z’ are Lebesgue-integrable in the 
neighbourhood of x’ = 0.: In general we are concerned with tie possibility 
that él,/ex’, I, may be unbounded at 2’ = 0; if they are bounded the 
analysis is considerably simpler. 

In the range of integration under consideration we may take % = 1 


and then - 


ig | 21 2) 
| [(z')b dz’ = I(0)log 7M . : y + O(1+2?). 
} a (n° oF 
—y 
7 
a es 
Hence : T(z’ )uos dz’ 21(0) | o(7). 
Ox . x x 
7 
Also 
€ 7] . 
lim | | Fol 2") = L(x’ ,z’) a4 dz'dx’ 
y—>0 J Cz OX 
0 —7 » * 
lim | | ( ub : he i. =) dz'dx' +-o | _ 
y>o J y Oz Ce fn 
0 —7 
are 
lim | (: A. +, 2 ) log{(a—a’)?+-y?\ da’ o( } 
y—0 , \Oz Cx M 
0 


where 
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We can split the range of integration [0, €| into the three intervals [0, 2/2], 





1/2, 3a/2], [3a/2,¢]. In the first interval we put y = 0 and we have to 
consider integrals of the type | ;- Let d(x) = | |1,| dz’ so that ¢(x) 
MH by MW 


; absolutely continuous and ¢(x) > 0 as «> 0. Then 


2 r/2 
da P dx 
] | Ft. 
~ i a 1 
0 0 
x/2 
» re 
( | L| dx’ 
x, . 
0 


In the interval [3a/2,¢«] we put y = 0 in the integrand and obtain an 


5 dx’ : 
integral of the type I, ——. We have 
32/2 
. i d. . x dx’ 
| ] t | L 3 t 
t a r av 
ar /2 3.r/2 
a Al 3X , x la’ 
; b(€) 2d, $4 | 4 ) oe. (1) 
€—a2 2 B ( x)* 


by an integration by parts. 
Now, given «, > 0, we can find X such that ¢(x’) < «, forO <a’ < X 
and we can choose x such that 2 = X®. Then, providing X < 3 and X <e, 


r X € 
- x dx P e, edz’ | h(a’ jar da’ 
D(X ) = ~4 r a . —j 
(2 he J (av’—z)® J (a —x/° 
3x /2 3z/2 x 
: (M—e,)at Ma 
1.€. a“€; - ; _ 
1—2z? €e—a 


where M is the upper bound of ¢(x’) in [X, e|. It now follows that as 2 > 0 
we can make the right-hand side of (1) as small as we please by a suitable 


choice of «,. Hence 
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Finally we consider 


3/2 , x 3x/2 ' 
v1 : dx el I . . dx’ 4- | aL ; : da 
, (~—x )*+y" ; " (x—2')?+y? } ~ (x—2' PF +y? 
r 3/2 
22756(3a/2)—d(a/2)} r F x—z2’ 
= ; z {P( n ) cal v/2)5 = | | { ed(22") c ; U : : dx'\ 
, x*-+-4y* 4 ox (x—x )*+Y* 
by an integration by parts. Now 
° a r—2’ P a 
ax h(x ) ee '\2 9 dx 9 9 (P(x) T 0 | I ); 
: ox (x—ax )*+-y* x*-+-4y* 


, 


since 4(x’)—d(x) o(1) as «+0, with 2’ in [a/2,2]. Similarly for the 
other integral. Hence 


37/2 
; 3 x—2’ ] 
lim I. he dx o{_). 
y0, (v—wa )°+-y" x 
r/2 
. \ 
é tims ite es I 
and thus lim | £ — log{(a—z2’)?+-y"! dz’ = o([-)}. 
y—0 , OX x 
0 
€ € 
— fat ak a , al oe ] 
Also lim | — log{(~x—wa’)?+ y?} da | — log|a—a2’| dx o( . 
yo J OZ J Oo v 
0 0 


Consequently we have shown that 


, 27(0) l 
B, i O () 


Thus in order that the boundary condition B, = 0 can be satisfied for 
small « we must have /(0) = 0, i.e. there is no excess electric current on (. 
Of the other possible excess line distributions mentioned, none provides 
a contribution from the near elements to B,, and so our conclusion remains 
valid even when the other singularities are present. 
3. The other line distributions 
We now suppose that the current perpendicular to C is non-zero so that 
a line of electric charge accumulates on C. The contribution of the near 
elements to Z,, is then 


€ 7) 


. 2 


lim 
y—0 - 


0 7) 


kpL,(a’,2)— 1 La',2') + La’, 2’)\ & 
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\ ea Cz CX - 





| dz'dx’ 4 
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|. b | LAO, 2 eb dz’. 
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[he dominant portion due to the electric charge is —2J/,(0,0)/x. The 
ntegral involving ¢l,/éx’+-cl./cz’ is of the type discussed in connexion 
with B, and is consequently o (1/2). The integral involving J, is clearly 
(1/~). Hence £, can vanish for small x only if (0,0) = 0, i.e. there is 


no line of charge on C and the current perpendicular to C vanishes at C. 


The remaining two line distributions, those of magnetic current and 
harge, provide no contribution to #,, #., or B, from the near elements 
and so they cannot be determined from the existing boundary conditions. 


The surface of a perfect conductor cannot support either magnetic charge 


current and so a reasonable and simple assumption is that the edge 
cannot support such distributions. This condition is equivalent to certain 
restrictions on the field. 


4. Behaviour of the field near the edge* 

We now introduce the cylindrical polar coordinates (p,¢,z),+ keeping 
the origin at O. The components £,, H., Bg are unaffected by the presence 
of magnetic current or charge on the edge. Hence the contribution of the 


near elements to By is 


eat T — 1, Ou 
- IT (x’,2 ) cos ¢+ L(2’, 2’) ; l ae dx’, (2) 
wy) . | C2 : op) 
u y 
where, in %, 1? = p?+a'2*—22’pcosd+2”. The integral involving J, has 


already been discussed in connexion with B,, for the cases ¢ = 0, 27 and 
shown to be o (1/p) as p > 0 in those cases. The first integral is known to 
be finite from potential theory. Hence By = o (1/p), and similarly we may 
show E, o (1 p). 

The contribution to #. from the near elements is 





€ T 
, ( pry c » ny) Ob a 
| [ [aeneew—1S Lez), L(x’, 2')\| de'de’. 
P - : (oa OZ J ez 
0 7 
It is not difficult to show that this expression is continuous across x = 0. 
In other words, FZ. is continuous across x = 0. Also since FE, = 0 for all x 
satisfying 0 < x < « it follows that FE, = 0 at x = 0, i.e. the component 


of electric intensity along C vanishes on C. Hence EF, = 0 (1) in the vicinity 
of the edge. 
Collecting these results together we have 
E o (1), E. 0 (1/p), Bg o (1/p). (3) 
* Note added in proof. For a recent discussion, from a different point of view, of the 


ehaviour of the field near the edge see J. Meixner, Ann. der Phys. 6 (1949), 2. 
' In this section ¢ has a different meaning from that in earlier sections. 
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When there are no line distributions of magnetic current and charge 


present we may easily show, by considering the contribution of the near 


elements, that 


B., = const.+-o (1), E o (1/p), B, 0 (1/p). (4) 


If magnetic current and charge were present there would be a term of 


O(log p) in B, and terms of O(1/p) in Ey, B, so that the absence of these 
singularities allows the behaviour in (4). Conversely, if we assume the 
behaviour in (4) the condition on E, ensures that no magnetic current is 
present and then the condition on B, requires that there shall be no 
magnetic charge. 


5. Uniqueness 

We can now show that the field is determined uniquely by (3), (4), the 
conditions on the surface, and the usual radiation conditions at infinity. 
Suppose that two solutions are possible. Let E, B represent the difference 
between them. The behaviour of E, B near the edge will be that indicated 
in the previous section. Let S,; be the surface of a sphere of radius R and 
centre in the neighbourhood of C, where R is taken large enough for S, 
to enclose S completely.’ Then 


$ (E\B*+E*,B).dS = 0, 


Si +S8e+Ss 


where E*, B* are the complex conjugates of E, B respectively. Now 
| EA B*.dS = 9 | (E, BYE, B)5 dbds 
S, om) 


o(1) 


as 6 > 0, from (3) and (4). Also 


| (EA B*+E*,B).dS = 0 

iS. 
since the tangential components of E vanish on S,. Hence allowing 5 > 0 
we have 

o (E, B*+E*,B).dS = 0, (5) 

Ss 
i.e. the integral of the Poynting vector over any closed surface enclosing S 
is zero. Consequently we may expect on physical grounds that E = 0, 
B= 0. We can show this mathematically by recalling that since E, B 
satisfy the radiation condition at infinity we can write, using spherical 
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olar coordinates R, 6, 4, 


EB — e-ikre NY (Ong 4 Gna gg Ong) 

— lr SP: hoa 

. (6) 
Bua eee | Pung, 4 Ome gy Ons i. | 

L \Ruvt tn * Res 


vhere i,, i,, i, are unit vectors in the directions R, 6, ¢ increasing respec- 


tively. The a’s and b’s depend only on @ and ¢, and the series are absolutely 


nd uniformly convergent for R > some constant (Ry, say). Also we know 


that i 
Ww , L — 
12 ] 913 13 k bis (/) 
nd that if a,, 0 and ay, 0 then E 0,B 0. From (5) and (6) we 
we for R A. 
f | ~~ a 9 = De ~ Ang — b*. 
124, Rn Z, Rr aa R" Ly R” 
¢ n=1 v=] n=1 n=1 
.a™ “. . at. + Oe. 
li Ane \” Sus Ans ‘, n2\ R2 dQ 0 
hat I™ Lee BR" a. R° 4 RY" 
} 1 ? 1 n=1 n=1 
vhere R? dQ dS. Hence, using (7), 
D (Ay2 47, +443 Ag) dQ 0. 
Since this is true for all S, such that R > R, we must have aj. = 0, a4, = 0 


for R R,, and hence E = 0, B= 0. Hence the field is determined 


niquely. 


6, A particular form for the current 
In this section we assume a special form for the current which satisfies 
ir previous requirements, and investigate the restrictions imposed by the 
undary conditions on the sheet. Suppose that near the edge we can 


xpand I, £, in the series 


iio .2 } vf f(z’) +a'f,(z’)+...}, 
Eig’. 2’) x'¥Sg, (z')+-a4'go(2z’)+...}, 
where a > 0, b lL, and f,, fo,..-. Jy, Jo... are continuous functions of z’ 


possessing continuous derivatives. We obtain from (2) that the contribu- 


tion of the near elements to B, is 


k2 er (, ( , 
“Ff, (O)cos@—a g(0 log(a’2— 2x pcosd } p”) dz’ +-(...), (3) 
LW Cp 
where f; (0) |df dz'|. _, and the dots in brackets represent terms whose 


form will be described presently 
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Now 


€ 


. 


w" log(w?— 2wp cos d+ p?) dw 


€ 
a+1 2? P st+1f p75 208 
am log(e?— 2ep cos d+ p?) a+1. 3 oe #, dw (8a) 
0 

by an integration by parts. Regard w as the complex variable | we and 
define w*+! by wt! w\atleta+l8 (0) < @ < 27). Let TI be the contour 
defined by the circles | w €, Ww 8° (0 <0 < 27, 8’ < p< e€) and the 
straight lines 6 = 0, 27 (8’ < |w| < e). Then 


~wetl(w cos d) dw ae = 
——= p = 2mrie4+7 o¢+1 cos{(a+1)(d—7)!. 
J w*—2wpcos¢d- p? 
- 

The integral over |w 5’ vanishes as 5’ + 0 and hence, when a is not 
an integer, 
€ 
"w"*l(w— pcos ¢) dw 

w?— 2wp cos d+ p? 

27 

mp"*1cos{(a+1)(d—z)}  e-Ha+D7—a+2 & (ei? p cos d)eia+2 do 


sin(a+1)z 2sin(a+-1)z J e®e?9—2epe" cos d+ p? 
0 


mp"*! cos{(a+1)(¢ 
sin(a--1)z 


7) j 


+ non-negative integral powers of p. 


We may proceed similarly for the term involving g,(0) in (8) when 6 is 
not an integer, and hence obtain for the contribution from the near elements 
2 
to Bs 
2ark*f; (0)p*+! cos{(a+-1)(4—7z)! 
: ; cosd 
ww(a-+1)sin(a+1)z 
: 2 b , 
2rk*q,(0)p” cos{(b+-1)(d—z)} 


: - ae F 9) 
tw sin(b-+- 1)z \ \ 


where the dots in brackets, like those in (8), indicate terms whose depen- 
dence on p is of the type p“*”"**, p?+"+! or p”, where n is zero or any 
positive integer. 

These bracketed terms of (9) include those of (8). The contribution of 
the far elements and of the incident field can also be included in the dots 
of (9) and so (9) represents By near the edge. 

In order that By = 0 when ¢ = 0, 27 the coefficients of the powers of p 
must vanish. Since neither a nor 6 is an integer. (9) shows that, if b 4 a+1. 


both a and 6b are odd half integers. If b a+1 we ean obtain no further 
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formation about a and 6 from (9) but the expression for E, is of the 
1 cosia(p—m)} 


. and 
S1In a7 


thus a and 6 are both odd half integers. 


If a is an integer, the integral on the right-hand side of (8@) can be 


valuated by elementary means to give a series of non-negative powers of p 


plus a term p“*!cos(a+1)¢log p. Making the appropriate modification to 


9) the condition on B, cannot be satisfied unless both a and 6 are integers 


such that b a+1. But then the condition on E, cannot be satisfied. 


9 


Hence all permissible values of a, b are given by the scheme 
a (2p+1)/2, b (2q¢+-1)/2, (10) 


vhere p, g are integers such that Pp 0,¢g 2 i 

It has been tacitly assumed in the above that neither /;(0) 0 nor 

0) = 0; this assumption is valid in the neighbourhood of some point 
f C unless one or other (or both) of the currents vanishes throughout a 
strip of width « with C as edge. The behaviour of the non-zero current 
vill be that found above. 

It is to be noted that, if in section 2 we had made the assumptions on 
he current used here, (9) and (10) would imply that there is no excess 
ne distribution of electric current on C; the form of £, then implies 


that there is no excess electric charge. 


7, Plane wave incident on a semi-infinite screen 

The analysis in the previous sections has been developed for the case 
vhere C is closed, but it is clear that it applies equally well if C extends 
to infinity as, for example, on a semi-infinite plane, providing that slight 
iodifications are made to S, and S83. 

We now consider the diffraction of a harmonic plane wave by a per 
fectly conducting semi-infinite plane. Wiegrefe (4) has indicated that a 
solution of this problem can be obtained by splitting the total field into 
two partial fields: (a2) that which has no component of the electric vector 
normal to the plane, and (b) that which has no magnetic component normal 
to the plane. It can be seen at once that the solution for (a) does not 
satisfy our boundary conditions at the edge since FE, = O(p-*). However, 
this case will not arise frequently since it can occur only when there is no 
electric charge on the plane. We shall obtain a solution satisfying our 
boundary conditions by means of integral equations. 

Let the edge of the plane be the z-axis and let the y-axis be normal to 
the plane so that the plane is defined by y = 0, 2 > 0. The direction of 
the incident wave normal is specified by the spherical polar angles 6,, ¢, 
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where 37 > 0, > 0,7 > ¢, > 0. The electric intensity of the plane wave | The V 


is given by (A, u, v)exp(—ax—By—yz), where equatio 
x = tksin 8, cos dy, B = iksin 8, sin 4,, y = tkcos6,. 

The constants A, », v are to be chosen suitably so that the wave is plane. : 

Let B,,, B,, be the discontinuities in the z- and z-components of B across | lim 

the conducting sheet and then, in view of the foregoing analysis, the as 


electric intensity is given by the equ 


E = (Mi+pj+vk)exp(—ax—By—yz)-4 It is 
part sO 
io ft PTs. ew oe 13) suc 
tae | eu B,.(x', 2’ )i—B,,(2’, z')k}+ Ie’) 
_ Eque 
c r , c , ’ , , , if the | 
f - B, Az’ ,2') ; B, A(x’, 2')\erad’ os | dz'dz', ; 
Cz Cx ‘ this ty] 
eaoe , : . , ; P It is, of 
where i, j, k are unit vectors in the 2, y, z directions respectively and in r bet ¢1 
: . ut the 
we put y 0 after the operations have been performed. B,,, B,, are , ; 
: ws a. : *, | of sim 
functions of 2’, 2’, but it is clear physically that their dependence on z ti 
, . equatle 
must be through e~”* and that only. Hence we may write a" 
unknov 
4rk? : ; 4k? : : -< 0: 
B,, = —J,(x')e-”*, BB, = ——— J {x')e-7”, it 
ig iw - iw” ofthe 1 
where J,, J, are functions of x’ only. We thus find, using Laplac 
r 29 , is 7 = and tr 
we-Y* dz mie-Y*H'?| e{(a—a’)?+-(y—y )?}#], 
that , 
7 ayo B ° "f 9 ! i c (2)T 6 1\9 2)4 , 
Ee”? = de-22-Bu ig | kT, yd, —*) NA cf(a—2' P+ yh] da’, 
, ei * Oa} Oa , 
0 
7 aye ° : CO c 9 "\o 1 ‘\o , : 3 
Bor = po-o2-B9—in | (yo, — 22) Halwa’ (yy VHH] de’, | ther 
0 , is cho: 
2 transf 
E. ev ,—an—By __ ; 27 ed, HN Sx V2 1 2b] da’ 
ye”? = ve Y—tr | (—y —F 1G [«{(a—a’)?+-y?}4] dx (e = ¢ 
0 , c¢> 


(ll) | o<r 
where H}?)(x) is the Hankel function of the second kind and 


k = (y2+k*)t = ksin6,. 
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The vanishing of Z, EZ, on the sheet now gives the following integral 
uations to determine J,, J: 


. 


(x? ye, “HY (x|\a—ax’|) dx oe Q, (12) 
. Cx 77t 
we Ae cJ.\ ¢ i de a. ‘ 
im 1.27. (J ) lH l\Ki(x~a—a )*-+ y?}* | da FT oe 0, (13) 
0 : CA CA } : 7 
the equations holding for a 
It is now convenient to assume that & has a small negative imaginary 
wt so that k k.—wk; (k, 0, k, 0). We search for solutions of (12), 
13) such that J,(x’) = 0 with a 0, J(a’) = 0 with 2’ < 0and such that 
x’) O(e-'), J.(x’) O(e-*’) for large positive 2’. 


Equations (12), (13) are simultaneous integral equations of the first kind 


fthe Wiener—Hopf type. The technique for solving single equations of 


this type is well known and a slight adaptation of it solves our equations. 


tis, of course, simple to transform (12) and (13) into two single equations, 
ut the following procedure is adopted since it will apply to the solution 
f simultaneous equations where the transformation is less simple. Each 


juation is extended to hold for all 2 by adding to the right-hand side an 


nknown function which is zero for 2 > 0 and equal to the integral for 
0: the term involving e-* is taken to be zero for x < 0. Then each 
fthe unknown functions is O(e'**) as a + —oo. We now take the bilateral 


x 


Laplace transform with respect to x, i.e. | e-*”, and assuming that the limit 


nd transforming integration can be inverted we obtain 


K2 g vs J, v ] 


: 2 . ys, (8), (14) 
(s?-+- x?) 2m s- 
ys (s? k?) j A ] “ 
é f Fy b,(8), (15) 
(s*+kK*) 2a 8 x 
vhere J. { J.(x)e** dx, Jf. | J(a)e-** dx and that branch of(s?+-«?)! 


0 0 
s chosen which reduces to « when s 0. Since k; 0, im(«) < 0 and the 
transform of the Hankel function is analytic in the strip im(x) < o@ < —im(kx) 
o+ir). Y,, # and 1/(s+a) are analytic bounded functions in 
re(a). w,(s), %(s) are unknown analytic bounded functions in 
re(x) (or —im(«)). Hence (14), (15) hold in the strip 


re(a) <o < re(a) 


since re(«) - im(x). 
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Eliminating Y, from (14), (15) we have 


i¢ yvs +A? 
271(s+-x) 


ysib,(8)+Kyh,(s), 





or ; 7 yvs+-AK? 
k*(s +4)! I2— 5 L x)(—a—ix)? 
ysihy(8)-+K%yfo(s)  ystAe® fo ol 1 |g 
(s—ix)! 2ri(s-+-a)\(s—ix)? (—a—ix)}!} | 
In (16) the left-hand side is analytic ino > —re(«), whereas the right-hand | 


side is analytic in o < re(a). Since the expressions have a strip in common 
it follows that both sides are equal to an integral function. Further, as 


s +0 with o > —re(a), the left-hand side is o(s!) since %,->0 as 
s —> 00; and as |s| > 0, with o < re(a), the right-hand side is 0 (s*) since 
b,(s) > 0, Yo(s) > 0 as |s| > 00. It follows, by the extension of Liouville’s 


theorem, that the integral function is a constant A. 


Hence | = A : f- ——— as cat — (17) 
k2(s+-ix)? § 2aik?(— a—ix)#(s+a)(s+ix)! 
and ysy(8) Facba(8) t a Beall é pia —} A. (18) 
(s—vk)? 2m1(s- x) l(s Ke)? (—a- -ix)*| 
But J,(2) > 0 as x > 0 and hence in (17) 
A4 a 0. 
271(—a—ix)! 

Therefore g, = (19) 


2mik?(— a—ix)*(s+a)(s- ix} 


Using the value of A in (18) we see that as |s| > 00, with o < re(a), ¢,(s) 
is o (1/s) or O(1/s) according as v is or is not zero. This was to be expected 
since ¥,(s) is the transform, over [—oo, 0], of 2,— FE, which is 0 (1) or O(1) 
as x > 0 according as v is or is not zero. 

We now return to equation (14) and observe that it can be written 
n° $.—ysJ, v(—a—ix)! v 


(s - ix), (8) +5 . 
“27r(s+-a 





m7 a {(s—ix)!—(—a—ix)}}. 
(s+7K)? 271(8+a) 

The argument proceeds as for (16) and both sides are equal to an integral 
function; using (19) and our information on the behaviour of 4,(s) at 
infinity, we find that the integral function is zero. Hence 


C ys(Ax?—vay) v(—a—ix)#(s+ix)! 


2rik?x?(—a—ix)*(s+a)(stix)? | 2rrik?(s-+ x) 


(20) 
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The application of the Mellin inversion formula to (19), (20) gives for 


J . 
Ak" a ° ola ik)u 
J Vr) ; / ie _ du. (21) 
2th? X ik) F u 
0 
d 
2aik? Y— 1K)? 
(16 f * o(ax—-tK)U 
, “)) KA Kf ayA—v(a?+ k?)te-™ — ~ 
u? 
‘ 0 
ani (22) 
nol 
as | In these formulae we may now take & to be real. In the particular case 
) as vhen 6; or, (21), (22) reduce to the usual formulae for two dimensions. 
nee The electric intensity at any point can be obtained by substituting (21). 
lle’s 22) in (11) but it is simpler to write (11) in the Mellin inversion form, e.g. 
I a he } expisx—i|y|(s?+-«?)!} ds, 
is s kK") 
(1% : 
vhere —re(a ry re(x). Substituting from (19), (20) 
S . 9 * 1 ° } 
EB \ l {| (vay Ak")(s ik)? yvs| X iK)? | 
: zi Jl e?(—a—ix)'(st-a) © x?(8-+-x)(s—ix)!] 
KEXP Sx d y (s?+ K2)h) ds. 
Dy sony | | (vay —Ak*)s 
MY, € i - - 
Jar . 1K=( v—ix)§(s-+ a)(s+ik) 
19 
yv{ x iK)2(s ix)? ; e 9\1 
: : | exp{sx—t y (s?+-«?)*} ds, 
1°) | K“(S—- ax) } 
ted 
M1 : v{ Y—tk) ‘ ; ‘ 
\ E.« ve ; ——— expjsx—tly |(s* K?)h) ds, 
27 § v)(S—tK)? 
where sen y l according as y 0. 
[It is to be noted that since the integrals for #,, E, are uniformly con- 
| vergent with respect to y we may put y = 0 in the integrand to obtain 
ral alues on the sheet. We see. therefore. that the inversion referred to 
at earlier is permissible. 
The behaviour of the electric intensity at large distances from the edge 
in be obtained by means of the saddle-point method of approximation. 
20) | Leta Reosd. y Rsind (0 < ¢ < 27), where R is large. Then the 
saddle point of sv—z\y|(s?+«?)? occurs at s ixcosd. In order to 
5092.12 
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deform the contour so as to pass through the saddle point it is necessary 
to pass over the pole s vif d, > d > 0 or 27 > o > 2x—4, but not 
otherwise. The region corresponding to ¢, > ¢ is the shadow and in this 
region the contribution from the pole removes the incident field. Corre- 
sponding to 27 > @ > 27—4, is the region of reflection and here the pole 
provides the reflected wave. The diffracted field is supplied by the saddle 
point, and for this portion of the field we obtain 


E. ev vay —Ak*?—K(Ak-+-ivy)cos 14! ikR 
4,€1% ~w : : SIN 5 
. (—a ix)!(a—ix cos d)K7! ‘s Rt 
er vy(a+ix)—K(Ak-+ivy)cos¢@ eWtkhR 
EB e?? mw — y \ Y) cos $4| sgn y - 


R} 


a*x(— a—ik)*(a—ix cos d) 
Soak v(—a—ix)! sin 46 e-kR 
E. Cr ~ = . z >1 7 
. im? (xa—ixcosd) PR 
except for values of ¢ near ¢, and 27—4,, for then the saddle point is in 
the vicinity of the pole s X. 
[ wish to express my appreciation of some helpful comments by Professo 
S. Goldstein. 
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NON-STEADY FLOWS UNDER ASYMPTOTIC 
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| SUMMARY 
In a previous paper (1) steady state solutions were obtained for the circulatory 
about a circular cylinder, igh the surface of which uniform suction is 
ied. 
In this paper approximate solutions are obtained in a closed form for the transient 
ecurring when the peripheral speed of the cylinder and/or the suction velocity 
ugh the surface is suddenly altered. The suction velocity is assumed to be 
ciently large for the asymptotic exponential velocity distribution to be closely 
iched in the initial and final steady states. 
[he velocity distributions in the transient state lie between those for the initial 
| final steady states and hence the circulation at infinity remains unchanged. 
5) effective duration of the transient state is inversely proportional to the square 
suction veloc 
solutions obtained for tl ular cylinder can be applied to the analogous 
lem of the non-steady flow past an infinite porous plate and for this case they 


1. Introduction 
) previous paper (1), solutions were obtained for the steady circulatory 
w about a circular cylinder, through the surface of which uniformly 
stributed suction was applied. In the general case, the cylinder was 
tated at a steady, arbitrary, angular velocity and there was an arbitrary 
rculation at infinity, but the suction had to be greater than a certain 
miting value for this circulation to exist. 
\ proposed method of generating this circulation was to rotate the 
linder with the suction off, thus setting up the well-known ‘free vortex’ 
w, for which the circulation is constant at all radii. Suction, greater 
han the limiting value, was then applied and the speed of rotation of the 
vlinder was then altered to the chosen value. 
\rguments were given to show that the circulation at infinity would 
main unchanged throughout this process and indefinitely. These were 
sed on considerations of the flux of vorticity by convection and diffusion 
ross a fixed circuit and depended on the hypothesis that, during the 
hange from one steady state to another, the vorticity was effectively zero 
itside a finite radius. 
It is desirable to establish the constancy of the circulation at infinity 
ore firmly, since this is of importance in connexion with the proposed 


[Quart. Journ. Mech. and Applied Math., Vol. 3, Pt. 4 (1950)] 
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use of the Thwaites flap (2) as an anti-gust device. In reference 1 it was 


shown how the Thwaites flap could be used as an alternative method of 


generating the circulation in the problems discussed in that paper. 

In this paper approximate solutions are obtained for the transient state 
occurring when the angular velocity of the cylinder and/or the rate of 
suction is suddenly changed. The suction is assumed to be always sufti- 
ciently large for the asymptotic exponential velocity distribution to exist 
in the initial and final states. These solutions are applicable to the infinite 
porous plate and are then exact. 

Since the appearance of this paper in its original form as A.R.C. 12040, 
B. Thwaites (3) has dealt with the same problem and has obtained a 


solution in series, whereas the present results are given in a closed form. 


2. Notation 
i radius 
q circumferential velocity component at radius 7 
a radius of permeable circular cylinder 
wy suction velocity normal to surface 
kinematic viscosity 


RR woalv 


t time 

Yq, peripheral velocity of cylinder 

Jo veference velocity K,, (27a) in case of circular cylinder 
q q do 

T vt/a? 

A roa 

y xr— |] (r—a)a 

C40) 4a 

K, circulation at r a 

K, circulation at 7 = a 


i: a constant 

N (R—2)/2 

qd, the part of 7 dependent on the time 
n distance from the surface 

Yo gniv 

T"  tq2iv. 


3. Theory 
3.1. The equation governing the motion 
As in reference 1, we assume the motion to be independent of @ the 
angular coordinate and a function of the time ¢ and the radius r. 
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The equation governing the circumferential velocity component q is (1) 


lL a o* l 
og t R 1) “1 1 (R 1) 4 (1) 
v ct or- r or ia 
ynere R =. (2) 
vy 


is the radius of the cylinder and wy is the suction velocity at the surface. 


Introduce non-dimensional variables g, x, 7’, such that 
Y 110 
ta 


which 7g IS a CONVeNnlent reterence velocity and is constant. Then 


1) becomes 


og 0-0 ? ] 7] S l - 
] I I ) og (R : Mg. (4) 
c Cx" a Cx } ee 
Assume a particular solution in the form 
q P k-T'4 R “st2). (5) 
Substituting in (4) we obtain 
d2 df 
ae co (k?a?— N?)f 0). (6) 
. R—2 = 
ynere NA = ) (1) 
4 
nd putting y, = && (8) 
6) reduces to the standard form of Bessel’s equation 
,d*f df ‘ ae 
yes vy (x27— N? 0. 9 
dx? Fe y (9) 
Hence a solution of (4) is 
q eKTy-K2f 4 J (ka)+ BY, (ka)}. (10) 


The steady state solution of reference 1 is also a solution of (4). This 


the general case was expressed as 
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7, is the peripheral speed of the cylinder and q, is given by 
K, 
0 . (13) 


2ra 


This we choose as our reference velocity. Hence 


" d = ( } (R-1) 
q a 1. | — L}a (14) 
do ( 
is a solution of (4). 


[f suffix (1) denotes the initial steady state at 7' 0 


l 
7 —1 »—(R,—1) 
q a | - 1): : (15) 
Vf 1 


and C, R refer to 0 < 7, then (14) will represent the final steady state. 
Equations (14) and (10) must be combined in such a manner that 


q=Ia | for all 7 6) 


q=h, T=0 
The building up of a solution to satisfy the boundary conditions and the 
end conditions appears to be complicated and we therefore proceed to 
simplify (4) by assuming R large. 


3.2. Solution for R large 

If & is large, we see from the graphs of reference 1 or equation (14) for 
the steady state solution that the rapid changes of velocity occur near the 
surface, x -> 1-0. Hence in (4) the term (R—1)q/2? is negligible compared 
with the remaining terms and (R+-1) may be replaced by R. Also x, where 
it occurs explicitly, may be replaced by unity. 


Put a l+y, (17) 
then (4) simplifies to es - ke Ro. (18) 
e7 ey? cy 

A particular solution of (18) is obtained as follows. Assume 
q = eH THR (y), (19) 
then (18) becomes a'f L qf 0. (20) 
dye * 
where x2 = k?—} R?, 21) 
Mquation (20) occurs in diffusion problems and a simple solution is 
f A’ cos a(y—B). (22) 


where A’ and £ are constants 


Putting 


F(B) 
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the solution of (18) is 


: F(p ’ 
q é tROT+I RY eos x(4— Pp) (23) 
Generalizing this, we may write 
7 | — ; a , 
q exp|— (1 R?7'+4Ry)} e-*T dy F(B)cosa(y—B)dB (24) 


0 
sa particular solution of (18), where F(8) is to be determined later from 
the end conditions via the Fourier integral theorem. 
The steady state solution, corresponding to (15), which is a solution of 
l : ° 
_ | | 1} Ry (25) 


Is is 


d for the initial steady state solution, 7’ 0, we take 


i=! (- 1); Rey, (26) 
C, / 


If C, R refer to 0 7’, then (25) is the final steady state, 7’ = «. The 
hange from ( to C represents a sudden change in peripheral speed (see 
\2)) and the change from #, to R represents a sudden change of suction 
elocity (see (2)). In the general case these can occur together. 

We therefore take as the solution of (18) the sum of the solutions (24) and 


2), 1.e. 


expi— (GAT } Ry); | ew dx | F(B)cos x(y—B) dB+- 
a ’ = 
L44 (o-1) Ry, (27) 
For 7 ©, this gives the final steady state (26). For 7 0,¢d=%.- 
Hence from (26) and (27) 


1 7 . 
0 t 


L} l + 1). Ke - thy | dx | F(B)cos x(¥/ B) dp 


1 "4 


0 


[——1) eas (=—1)ea™ {da [ F(B)cos «(y—B) dB. 


(28) 
Hence by the Fourier integral theorem 


F(B) - 1) Ry) RB (7.1) dS (29) 


where, so far, F(B) is only defined for 0 - B < 0. 
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The boundary condition at the surface requires ¢ = q,, y = 0 for all 
T > 0, i.e. Gg = 1/C. Hence from (27) 

F(B)cos xB dp Q, (30) 

This is achieved, since cos xf is even in 8, by making F(8) odd in 8. Hence 


for 
0 <6 1) 


F(—B) F(p) 
where F(f) is given by (24) for 0 < B 0. 
Hence (27), together with (29) and (31), is the required solution of (18), 


satisfying the boundary conditions and the end conditions. 
Writing 


(31) 


. I - o . 
dr = —exp{—(gFPT + $Ry)} | e-8F do | F(B)cosaxa(y—B)dB (3: 


0 x 


and assuming that the order of integration can be changed, we integrate 
with respect to «, using the result 


e* cos 2hz dz e-l* (33) 


e-GQRT+tRy) ¢ 


to obtain Gr - . F(B)e (y—BP AT dp, (34) 
2, (7T) 


which, since F(8) is odd in B, can be written 
ek RT +4 Ry) 


ae p (y—B/47'__ p(y +BP/AT’) , 
1 let | F(B)(¢ ) dp 


0 


V7 


URT4yT~ PP 
aaa abil , ddl 1) (Ry R)B 
2i(rT) AC, 


) 


on substituting for /(8) from (29). 
Multiplying out the brackets and putting 
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35>) gives, on rearrangement 


eI / 1 : 
7] ef ( I ap el er D) 1p 
11 a ite Ip, | iby 
( 0 
, " 
(- l . e~ 1 dB, eb | Lf By ia| ( v ) 
0 4 
here 
A 1(RT*—yT-*) }( R?T)-*( R?T — Ry) 
B 1( RT! +-yT-?) i( R°T)-*( R?T + Ry) 
» 
D 1f(2R,— R)T*?+yT-} 1(R?T') (2 1) R274 Ry| >, (38) 
: \\ R a 
»p 
E LS(2R,— R)T?—yT-} 1(R*T') f(A 1) R27 Ry\ 
2 2 \\ R | 
The integrals in (37) can be expressed as error functions. Thus 
f a | | ry en D? : 
TF - {Fa 1) | (l1—erf H)—e”’'(1 erf D)| 
a Ife" erf A)—e”"(1—erf B)) (39) 
_ ‘ l ' 
nd so from (27) q qy ‘ G i) ~~}. (40) 


Check. That this is a solution of (18) is easily checked, since each of the 


terms for 7; in (39) satisfies the equation. 


Putting y = 0, makes q, vanish and so from (40) 
i l 
q=7 ie. I=d as required. 
Putting ¥ © makes g, vanish and again from (40) 
qd l, ie. G¢ = as required. 
Putting 7’ © we are left with the final steady state solution 
l > 
q ] + | ] é Ry, 
ra 
Putting 7’ — 0, after multiplying through by e-”", in (39), gives 
77...) 1) %e-Riy _ 1 }2¢ Ry\ 
lv) T-0 = =f “ a 
2\\C, ( 


ind substituting in (40) we obtain 


ws required. 








442 





J. H. PRESTON 


The velocity distribution given by (38), (39), and (40) is seen to be a 
function of Ci, C, R,/R, R?T', and Ry. The C, and C refer to the initial 
and final peripheral velocities of the boundary: R,/R relates the initial 
and final suction velocities, for, from (2), 


2 @ . 
R, (Wo),@ v (Wo), 4] 
R » wa uw rae) 
t t 0 0 
t?7' is, from (2) and (3), given by 
wea vt ut 
R?T > » 40 (42 
vy @ v 
[f » denotes the distance measured from the surface. then 
r—a n ; 
y= ; (43) 
ai a 
w,an wh 
so that Ry : —, (44) 
V a V 


Hence, we see that R enters the problem only in ensuring that the 
exponential velocity distribution is attained. When this is so, the radius a 
of the cylinder is of no significance and the essential non-dimensional 
parameters, describing the velocity distribution ¢, are 

0 o 0) wet Wyn 


" ; ’ and 
uw Vv Vv 


3.3. Particular cases 


R, R: suction kept constant. 
(a) C = |: peripheral speed to generate appropriate circulation, i.e. 
( 
da = U at start, T = 0. 
C 4%: peripheral speed zero, i.e. g, = 0. T > 0. 
Then _ (RT) 
A= —— (7 Ry) 
nf (45) 
pain 
B = D = ‘=! (RT + Ry), 
and from (39) and (40) 
d 1fe-ky(] —erf A)—(1—erf B)!+(1—e-¥»)., (46) 
(py | R, ~ R: sudden change of suction flow at T = 0. 
) 
(, = C = a: eylinder at rest at all 7. 


Then A, B, D, and £ remain as given by (38), and from (39) and (40) 
eB , : : . ; a sos 
q - e+ (1—erf A)—e”*(1—erf B)+-e”"(1 —erf D)—e¥*(1—erf E)! 


+(1—e-*¥), (47) 
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If R, ©: infinite suction, then (47) reduces to (46) with A and B given 
yv (45). 
4, The infinite flat plate with uniform suction 

Keeping to the present notation, let wy be the suction velocity normal 
to the plate, g be the velocity parallel to the plate at a distance n from the 
surface, and g, the value of g for n — oc. Here gq is assumed to be a function 
fnand?t only. Then the equation of motion is 


9 
CY Cc cf 

V : 7 Wo g (48) 
ol on- on 


Since there is no representative length associated with the plate, we 


itroduce non-dimensional variables defined by 


q W1 
, Vv 
n Y \ 
Yo { (49) 
t (id . 
10 
nd (48) becomes Ta w Cg 
or “7 
4 c Ui :; P a. (50) 
cl cy ~ Yo] CY 


This is of the same form as (18), the boundary conditions are the same, 
ind so the solution—equations (38), (39), and (40)—will hold if we replace 
R by (wy Go), T’ by T’, and y by y’. The non-dimensional parameters given 
n equations (41), (42), and (44) still stand. 

Since equation (48) is exact and no approximations have been made, the 
solution found for the cylinder and applied to the plate is exact. This 
rises from the fact that the exact steady state solutions take the exponen- 


tial form in the case of the infinite plate. 


5. Example 

As a simple illustration of the change of velocity distribution with time. 
we take g to be given by (46), section 3.3. This velocity distribution arises 
in two different problems. In the first, the suction is kept constant, but 
the angular velocity of the cylinder is suddenly reduced to zero from the 
value necessary to generate the circulation at infinity. In the second, the 
cylinder is at rest with a finite circulation at infinity and the suction is 
suddenly reduced from an infinite value to a finite value. Corresponding 
flows exist for the infinite flat plate; ¢ = q/q» is then a function of w%t/v 
and wynjv. In Fig. 1 ¢ is plotted as a function of wyn/v for a few chosen 
values of wt/v. When wt) 9-0, g is everywhere within 0-002 of the 
steady state value ¢-> oo. It will be noted that the velocity distributions 


lie between the initial and final velocity distributions. 
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In the case of the cylinder R = w,a'v > 100 for the exponential form 
to be attained (see reference 1). To take a numerical example, with 
wet/v = 9-0, if a L ft., v = 1-59 10—+ ft.2/see., then Ww, > t ft./sec. for 
R > 100, For wy = } ft./sec. t = 9-Ov/w2 = 0-09 sec. 


——_——_—_—_——, 


0-3- 





0:8 












































0-6 pine ‘mae ween 
0:5 iia 
ga / 
“Yo 
0-4 poe 
.% = Ke/2z | | 
| Wo =Suction velocity normal to surface | 
0:34 > n = Distance from surface —- _ 
, et Poy ¢ 4 
| | | | 
0:2 +— , — —__—__+—— ————$_____}____} _——— 
| 
Q it a ae ee T Se aaa t nee 
on 
a eee = 1 
0 1 2 3 4 5 a 8 9 


Wp 7/V 
Fic. 1. Velocity distribution q¢/qg v. wan/v at various times ft, (46). for (a) 
rotating cylinder with steady circulatory flow : suction applied, then cylinder 
suddenly brought to rest ; or (6) cylinder at rest with steady circulatory 
flow and infinite suction, then suction suddenly reduced to finite value. 


In the case of the infinite plate, if we take a free stream velocity 
qo = 200 ft./sec., then the times for the final steady state to be effectively 
reached for various values of w/q. are as follows: 


w 


: 10-3 w'(ft. /sec.) 0:2 t(sec.) = 0-036 
do 


10-4 0-02 == 3-6 


10-° 0-002 = S00. 














NON-ST 


6, Con 
Prov 
the tra 
eylinde 
The 
state Si 
that tl 
The 
referel 
could 
referel 
edge, 1 
into tl 
incide 
scarce 
The 
cable 
The 
the s¢ 


short 





for m 


with 











VON-STEADY FLOWS UNDER ASYMPTOTIC SUCTION CONDITIONS 445 


6, Conclusions 
Provided the suction is sufficiently large, solutions can be obtained for 
the transient state occurring when the angular velocity of the circular 


vlinder and or the rate of suction is suddenly changed. 





The velocity distributions lie between those of the initial and final steady 
state solutions and the circulation at infinity is unchanged. It seems likely 

it this will be so for smaller suctions, provided R (= wya_v) :. 

The proposed method of generating the steady flows described in 
eference 1 should therefore work, and it is highly probable that they 
uld be generated by means of the Thwaites flap as also described in 
eference 1. If this is so, then in the case of a wing with rounded trailing 
ize, the Thwaites flap, having generated a circulation, can be withdrawn 
nto the wing and the circulation will remain independent of velocity and 
ncidence changes. The lift will depend only on the forward speed and will 
urcely be affected by gusts, as was pointed out by Thwaites (2). 

The solutions obtained for the circular cylinder with R large are appli- 
ible to the infinite flat plate for all suction velocities and are exact. 

The effective duration of the transient state is inversely proportional to 


he square of the suction velocity and for normal rates of suction is quite 


L owt 
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FLUID MOTIONS WHOSE KINEMATICS ARE 
INDEPENDENT OF THE COMPRESSIBILITY 
OF THE FLUID 


By D. H. PARSONS (Exeter College, Oxford) 
[Received 29 November 1949] 


SUMMARY 

The following question is studied: What steady motions of a compressible fluid 
may be reproduced exactly by the motion of an incompressible fluid ? More precisely, 
what vector functions of position exist, representing simultaneously the velocity in 
a possible motion of a compressible fluid, and the velocity in a possible motion of an 
incompressible fluid? A general account is given of the theory of motions of this 
type. Two-dimensional motion is then studied in detail, and it is shown that in order 
that a two-dimensional fluid motion may be of the required type, it is necessary and 
sufficient that the streamlines be either concentric circles or parallel straight lines. 
A very special case of three-dimensional motion, in which the streamlines are parallel 
straight lines, is then mentioned. Finally a general discussion is given of the type 
of results obtained, and it is shown how the restrictions encountered are interpreted 
by means of the theory of partial differential equations. 


1. General theory 

We consider the steady motions of two non-viscous fluids; firstly a com- 
pressible fluid, whose pressure is related to its density by the adiabatic or 
other prescribed law; secondly a uniform incompressible fluid. In either 
case let p, p, u denote respectively the pressure, density, and velocity 
vector at a point whose position vector is r. Then the object of this in- 
vestigation is to determine in what types of motion, if any, the same vector 
function of position u(r) can represent the velocity in both cases. Thus 
the problem is to find a single vector function u(r), such that the equation 
of continuity and the dynamical equations of motion for either fluid may 
be satisfied by a suitable choice of the (scalar) pressure function p(r) in 
each case. For the sake of brevity we will henceforth refer to the com- 
pressible fluid problem as case (a), and to the incompressible fluid problem 
as case (b). 


In either case the equation of continuity is 
div(pu) = 9, 
which in case (a) may be written 
pdiv u+u.grad p Q), (1) 


and in case (b) becomes div u 0. (2) 


[Quart. Journ. Mech. and Applied Math., Vol. 3, Pt. 4 (1950)] 
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FLUID MOTIONS 





Since U must satisfy both (1) and (2), we have in case (a) 
u.grad p 0, (3) 


showing that the density is constant along any streamline. Now the 





| dynamical equations of motion, in either case, may be written vectorially 


| S 


grad) }u? | (] p) dp, u A curlu. (4) 


Forming the scalar product of this equation with u we obtain in both cases 
u.grad(}u*)+-(1/p)u.grad p = 0, (5) 

hich in case (2) may be written 
u.grad(}u*)+ (1 p)(dp/dp)u.grad p 0. (6) 

[hus from (3) and (6) we see that 
u.grad(}u*) Q), (7) 
owing that the velocity is constant in magnitude along any streamline. 
Next, taking the curl of (4), we obtain in both cases the well-known 

ecessary condition for a steady fluid motion, 

curliu ( curl u! 0, (3) 
nd we will now see that the conditions (2), (7), and (8) are sufficient for 
the solution of the problem. For suppose that u(r) has been found, 
satisfying these three equations. Then (8) is the sufficient condition for 


the existence of a scalar function of position, x(T), such that 


uA curlu grad y; 


nd hence we may solve (4) for | (1/p) dp in terms of the position vector r. 

Thus ; _— , 
| (l/p) dp x }u°—+constant. (9) 

In case (a), | (1/p) dpisa known function of p, whose derivative (1/p) dp/dp 


snot zero, and hence (9) gives p, and therefore p, in terms of r. While 
nease (b). | (1/p)dp = p p)+constant, py, being the fixed density, so that 
nee again from (9) we may find p in terms of r. In each case, (4) being 
now satisfied by the appropriate choice of the pressure function, we obtain 
5) as before, and hence in case (a), (6). From (6) and (7), which is by 
ypothesis satisfied, we deduce (3), in case (a); and (2) being also satisfied, 

1) then follows. Thus all the conditions in both cases are fulfilled. 

There are several applications of the foregoing theory in which large 
simplifications of the conditions obtained exist. In motion in three dimen 
sions, with symmetry about an axis, (2) is automatically satisfied by the 
se of Stokes’ stream function Y. If r be the radial distance from the 
ixis, and ¢ the vorticity, (8) gives the well-known result that ¢/r is constant 


long any streamline; and hence that ¢ r is a function of ‘. The condition 
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(7) then expresses the fact that the magnitude of the velocity is a function 
of ’; from these last two conditions, therefore, two partial differential 
equations, one of the second order and one of the first order, are obtained 
for the determination of ’. 

In two-dimensional motion the simplifications are much greater; and 
this problem, and also a special case of motion in three dimensions, will 
now be treated in detail. 


2. Two-dimensional motion 

Let u, v be the components of velocity parallel to the rectangular axes 
of x, y respectively, and let ¢ be the vorticity, so that 2¢ = (@1/éa)—(éu/éy) 
Then (2) is satisfied by the use of the stream function J, so that 

u c us cy; ‘i 4 us Cz. 2¢ V-u: 
while (7), which exhibits the constancy of the velocity along a streamline, 
shows that the velocity must be a function of 4. Thus 
(Cus /Ca)?4 (Cys cy)? function of us. (10 

In this case (8) becomes a single scalar equation, 

uclicx+v el/ey+Cleu/ex+-ov/cy!} 0, 
or using (2) ' wolj/ex- vecicy 0, (11) 
which expresses the standard condition that the vorticity is constant along 
any streamline. 

Now (10) is the well-known necessary and sufficient condition that the 
curves ¢(x, y) = constant, in this case the streamlines, shall be a family 
of ‘parallel’ curves; that is, curves having a common system of normals, and 
such that the distance along a common normal between any two chosen 
curves of the parallel system is constant. Let us then suppose that a 
family of parallel curves has been chosen. We now investigate the vorticity 
of motions having this set of curves as streamlines. Let [ denote a certain 
streamline; and let @ be the angle between the x-axis and a tangent to I’: 
and let Y be the (constant) normal distance between I and a selected 
reference curve of the parallel system. Then the equation to the tangent is 

«sin 6—ycosé— F(@)—Y = 0, (12) 
in which F(@) is a function depending on the form of the parallel curves. 
Then except when @ is fixed, in which case the streamlines are parallel 
straight lines, [ is the envelope of (12) as @ varies, and is thus given by 


» ‘ 
(12) and x cos 0-+-y sin @— F’(0) = 0; (13) 


which is clearly the equation to one of the common normals to the parallel 
curves. The coordinates of the point of contact of the tangent to I, i.e. the 
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tersection of the tangent and normal, are given by solving (12) and (13) 


y and y. Alternatively we may regard (12) and (13) as expressing Y 


|@in terms of x and y: hence differentiating and using the equations 


emselves we obtain 


(oY /ex)? (oc) cy)= g 
e?Y /x*) + (e*Y /oy? 1 {Y+ F(6)+-F'(6)j. J 


constant, there must be 


(14) 


t since the streamlines are given by (x, y) 
nique correspondence between % and Y; so that % must be a function 
Y. such that di /dY 0. Thus 
—= Vv (dys dY)V?Y 


(d*b/d¥* 


(d*b/dY*){(eY /ex)* 
(dis dV) SY 


erefore in order that the condition of the constancy of ¢ along a stream- 


(¢ r cy)*} 
rv (14 2¢ F(0)+- F’(@)}. (15) 


may be satisfied, that is, in order that (15) may be independent of 6, 


s necessary and sufficient that 


F"(0)+-F(0) = ¢, (16) 
ere Cc is a constant. The general integral of (16) is 
F(@ g sin 6—f cos é+ce, a 
: mii (17) 
F'(6 q cos 6+-f sin 6, 


ere g and f are also constants. Then substituting from (17) in (12) and 


3) we have for the tangent and normal respectively 


(a —g)sin 6—(y—f )cos 6 Y +c, (12’) 

(a—g)cos 0+-(y—f )sin@ = 0. (13’) 
us squaring and adding we obtain 

(x—g)?-+(y—f)? = (Y+e)2, (18) 


wing that the streamlines are concentric circles, the centre being the 
rhitrary point (g,/ ) 

In the special case when @ is fixed and the streamlines are parallel 
‘ight lines, given by (12) 


(CY /éa 


we obtain directly from (12) 
| (c?Y /ex?)4 


ds /dY2, (15’) 


(¢ Y /ey)? (c?¥ cy?) 0, 


| thus 2 


wing that ¢ is indeed constant along a streamline. 

Thus we see that the necessary and sufficient condition for a two- 
iensional motion of the type required is that the streamlines shall be 
ther concentric circles or parallel straight lines. Important special cases 
te the irrotational motion generated by a line-vortex. and the motion 


roduced by a vortex of circular section. 


99 ] 
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3. A special three-dimensional motion exce] 
We now examine the situation which arises in three dimensions when lepe 


the velocity vector is fixed in direction, which, of course, implies that the | acco 


streamlines are parallel straight lines. Let i, j, k be three orthogonal unit Sil 
vectors, parallel respectively to the rectangular axes of x, y, z; and let the | may 
velocity be always parallel to the axis of a, so that ecti 


. fund 
u U(x, Yy,2)1. “ 
de 
Then (2) gives 6u/éx = 0, showing that the velocity is a function of y and | jocc 


> only; and (7) is then also satisfied. For the remaining condition we have 

curlu = (éu,/éz)j—(éu/ey)k, 
u A curlu (u Cu/ey)j+(u éu/ez)k 
grad(}u’*), 
and thus (8) is always satisfied. Thus, for example, any steady flow of a 
gas in a straight pipe, in which the velocity varies over the cross-section 
in any manner, may be kinematically reproduced by the flow of a liquid 
in the pipe. 


4. General conclusions 

It seems at first sight rather surprising that the results of the investiga 
tion into the two-dimensional case should be so restricted. We have in fact 
seen that the possible streamline systems are obtained by equating to a 
constant the simple function on the left-hand side of (18), which contains 
only two arbitrary constants. Now the condition (8) applies to all stead) 
fluid motions and (2) applies to all motions of an incompressible fluid 
Furthermore, although we have added the condition (7), the motion has 
not been restricted to being irrotational. Thus the results might have been 
expected to be of some generality. 

The explanation lies in the fact that while the condition for irrotational 
motion, curl u = 0, involves (8) as a direct consequence, the condition (7 
does not. If, with the notation of section 2, we substitute for w and v in 
terms of ys in (7), we obtain a partial differential equation of the second 
order for %, of which (10) is an intermediate integral. Again. if we substi 
tute for uw, v, and ¢ in terms of % in (11), which is equivalent to (8), we 
obtain a partial differential equation of the third order, of which the 
condition of the constancy of ¢ along a streamline is an intermediate 
integral. Thus fundamentally the two-dimensional problem is that o! 
finding the common solutions to two partial differential equations, one 
of the second and one of the third order; and it is a well-known result in 


the theory of differential equations with two independent variables, that 
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xcept when special conditions are satisfied, the solutions to this problem 


when | depend upon not more than six arbitrary constants: this is entirely in 
at the cordance with the result obtained directly in section 2. 

iL unit Similar remarks apply to the case of motion with axial symmetry. It 
et th ay readily be seen that the conditions mentioned briefly at the end of 


ection 1 are intermediate integrals of the equations which express the 
ndamental conditions. These are once again of the second and third 
rders, so that just as in the two-dimensional case we should expect the 
'= sible solutions to be of an extremely simple and restricted nature. 


estlg 
in tact 
g to 
ntains 
stead 
fluid 
on ha 


e bee! 


itiona 


ion 











THE CHARACTERISTICS OF SYSTEMS WHICH 
ARE NEARLY IN A STATE OF NEUTRAL 
STATIC STABILITY 


By W. J. DUNCAN (College of Aeronautics, Cranfield, Beds.) 


| Received 14 Mareh 1950] 


SUMMARY 

It is shown that the rate of subsidence or divergence A of a system which is near 
a state of neutral static stability can easily be calculated trom a knowledge of thi 
mode of displacement in the neutral state and this mode is found by solving a set 
of linear algebraic equations. The first-order correction to the mode ean also bi 
found and in important cases this can be made the basis for calculating a second 
approximation to A; if necessary a further correction to the mode can now be found 
and from this a still more accurate root can be caleulated. The method can bi 


extended to continuous systems having infinitely many degrees of freedom. 


1. Introduction 
SysTEMs which are nearly in a state of neutral static stability are of con 
siderable interest and technical importance. Exactly neutral systems are 
specially simple inasmuch as the neutral mode of displacement can be 
found merely by solving a set of simultaneous linear algebraic equations. 
It is shown in this paper that the small rate of divergence or subsidence ) 
of a nearly neutral system can be very easily calculated to the first orde: 
of small quantities from a knowledge of the neutral mode; the first-order 
corrections to this mode can readily be found also. When the matrices of 
the inertias, dampings, and stiffnesses are symmetric the mode correct to 
the first order can be made the basis of a caleulation of A correct to the 
second order. For a continuous system with infinitely many freedoms 
the first-order approximation to A can be derived from the values of the 
potential energy and rate of dissipation of energy in an imaginary divergent 
motion proportional to the neutral mode and having A equal to unity. 
2. Exactly neutral systems 

The matrices of the inertias, damping coefficients, and stiffnesses are 
as usual, denoted by A, B, and C respectively. These are square matrices 
of order n, where x is the number of degrees of freedom of the system. 
The static displacements ¢ corresponding to static loads Q are given by 

C'q 0 Za 

or q C-"0 
If the determinant of C becomes very small, the elements of C~! become 


[Quart. Journ. Mech. and Applied Math., Vol. 3, Pt. 4 (1950) ] 
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ery large and the dis}'a‘«) ents g corresponding to given loads Q also 
come very large. Whe vanishes the system is in a state of neutral 


static stability and then ‘nite displacements can occur in the absence of 


pplied loads Since the constant term in the determinantal equation 
fthe system is equal to C', the equation has the root A 0 when |C 
wishes. Thus the velocities ¢ and accelerations ¢g are zero in the mode 
f displacement corresponding to the root A 0. There are, however, 

1) other modes which, in general, correspond to non-vanishing roots A, 
nd in these modes the velocities and accelerations are not, in general, zero. 
ses of multiple zero roots can arise but are rare. When a displacement 


curs in the neutral mode the forces are always balanced. Hence the total 


ork done in a neutral displacement is zero. As an example, when the 
ngle of incidence of an aircraft varies slowly when the centre of mass is 


t the longitudinal neutral point the resultant aerodynamic force is con 


stantly equal to the weight of the aircraft. 


> 


3. Nearly neutral systems 

Suppose next that the system is near but not in the state of neutral 
static stability. There are then some simple methods for finding the 
aracteristics of the nearly neutral mode. 

a) The matrices A, B, and C are all symmetric 

In this instance we can apply an extension of Rayleigh’s Principle. 
Suppose that in a datum state of the svstem there is a characteristic root A 
nd a corresponding modal column g. Then A must satisfy exactly the 
jwadratic equation 

A*(q' Aq) +A(q' Bq)+q'Cq 0, (3.1) 
here the coefficients are quadratic forms in the elements of g. This 
juation is obtained from the dynamical equation in matrix form 
\*Aqg+ABq+Cq 0 

v premultiplication by qg’, the transposed of q, i.e. the row matrix having 
e same elements as the column q: the equation is true whether A, B, 
nd C are symmetric or not. Now let A, B, and C be symmetric, and let 
he elements of g depart from their true values by small quantities of the 
rst order. Then according to the principle the value of the root A obtained 
rom the quadratic equation (3.1) will only be in error by a small quantity 
{the second order (1). Accordingly, let A, B, and C refer to the near 
eutral state and substitute for g the neutral modal column gq» which is 
sily caleulated. The root will be a small quantity of the first order and 


will be of the same order as the terms neglected and therefore can be 


We ppose tl n started by an extremely small impulse. 
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discarded. Hence we derive the approximation 


1 = — 10% (3.2 
4o Bq 
which is easily computed. 
Let C reduce to C, in the exactly neutral state so that 
Cod 0, (3.3 
and let also C = G,+8C, (3.4 


where the elements of 6C' are small quantities of the first order and (, is 

symmetric. Then by (3.3) and (3.4) 

Jo CAM 1 OC JW (3.5 

r GW OC JW 
Yo Bq , 

which is correct to the first order of small quantities. From knowledge of 


and (3.2) becomes (3.6 


this approximation to A we can calculate the first-order correction to the 
modal column qd» merely by solving a set of simultaneous linear equations 
(see equation (3.13)). 

The approximation to A can now be improved by use of the corrected 
modal column and the equation (3.1). Since A is a small quantity of the 
first order, equation (3.1) will be correct to the second order if q' Bq is 
correct to the first order and q’Cq to the second order. Now suppose that 
all the elements of 5C are proportional to p, which is thus of the first order. 


and let Na nC. (3.7 


Suppose also that to the second order 
q dot PN Pq. (3.8) 
where q, is now known but q, is not. Then to the second order of small 
quantities 
GCF = PQ ONG I) FPG oF HANFHOADELNG I) (3.9 
in view of (3.3). But equation (3.3) gives by transposition 
qo ( 0 0 (3.10 
since () is symmetric. Hence (3.9) becomes 
GCF = PWAG+ PVD AG+N ONAN G Wo): (3.11 
Now this does not contain q,, and we conclude that equation (3.1) will be 
correct to the second order of small quantities when we substitute for 4 
the approximate modal column which is correct to the first order. Hence 
the small root of (3.1) will now be correct to the second order. 
The improved approximation to A obtained from the quadratic equation 
(3.1) can be used to obtain an improved approximation to the correspond- 
ing mode by use of the dynamical equations. Then this recalculated mode 
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iy, if desired, be used to obtain a still better approximation to A by 
eans of (3.1). The whole process can be repeated as often as desired. 
b) The matrices A, B, and C are general, exce pt that C is nearly singula) 
The dynamical equations in matrix form yield 

(A°7A +AB+C,+6C)(qg+5q) = 9, (3.12) 


)» 9 


hile (3.3) is again satisfied. If we now neglect second-order terms we 


btain ' » 34 

Abd, Cy og oc do: (3.13) 
is represents » linear scalar equations and there are apparently (n-+-1) 
nknowns, namely A and the x elements of 6g. However, the modal column 


sarbitrary to a scalar multiplier and we may accordingly assign one of the 


elements of dg, say make one element zero. We have now n unknowns, 


hich can be caleulated 
We shall now show that when (, is symmetric this procedure leads 
exactly to the value of A given by (3.2) or (3.6). Premultiply (3.13) by qo- 


We obti Set ' 
" — Ado Bq dol 9 09-44 OU qo) 
% oC do 


(3.10) and this is identical with (3.6). 


4, A simple numerical example 


Let the inertia matrix be 








(.2 4 
1 as aS (4.1 
La oS eg 
the damping matrix 
2. 29 
b _ a os (4.2) 
> ee 
id the stiffness matrix be 
C = 480, (4.3) 
vhere 
a a 
( 3 6 | (4.4) 
a 2. 
ind is singular, while 
Oo VO O 
d¢ p| 0 0 Of, (4.5) 
00 14 








vhere p is supposed to be small. The matrices are all symmetric so both 


methods described above are applicable. To facilitate comparisons with 
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the approximations which we shall obtain we give the expansion of the 


exact determinantal equation for A, namely 


9 


WA+AB+C S3A8 +- 96A? + (16 p+ 144)A44 

(l4p+86)A3 + (21 p+-33)A2+-(8p+-6)A+p = 0. (4.6) 
When p is zero this has a zero root in accordance with () being singular. 
The modal column q, corresponding to the zero root is obtained at once 


from the linear equation ' 
(odo = 9 


and is found to be 
») 


Yo l : (4.7) 

] 
where the last element has been assigned the value unity. When p is given 
the value 0-1 we find from (4.6) the true value of the numerically small 


root and the corresponding modal column to be 


A 0-O1597, (4.8 
2°22030 
q 1-13840 (4.9) 
1-0 


to five places of decimals. 
The first approximation to the small root as given by (3.6) is 
A Lp 0-01667. (4.10) 
Equation (3.13) yields the same value of A and when we assign the zero 
value to the final element in 6g we obtain the approximate modal column 
2-21667 
1-13333 ; (4.11) 
1-0 
On substitution of this value of ¢ in equation (3.1) we obtain the quadratic 
equatio Se ee 
jeawon 22-12778A7+-7-20500A +-0-10944, 0, (4.12) 


of which the smaller root is 


A 0-O1597,. (4.13) 


This differs from the true root only by about 1 part in 5,000, The error 


in the modal column is of the order of } per cent. and the departure from 
the neutral modal column is about 10 per cent. 


For comparison, the calculation has been repeated with the inertia 


matrix reduced to one-quarter of that given in (4.1), i.e. with all the 


inertia coefficients reduced to one-quarter of their former values. [he 


first-order approximations to A and to the mode remain as in (4.10) and 


4.11 
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appr 
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of the $11) respectively. However, the quadratic equation (3.1) giving the 
econd approximation to A is obtained from (4.12) by reducing the 
efficient of A* to one-quarter of the value in (4.12) and the second 


pproximation 1s d 0-01537,. (4.14) 


ul; Yow the ‘exact’ root as obtained from the determinantal equation is 
C one A 0-O1537,. (4.15) 


so the error is only about 1 part in 15,000, or say one-third of the former 
rror. The ‘exact’ mode is 

2-22377 
+7 q 113745 |, (4.16) 
1-0 


give! 0 the first approximation to the mode is in error by roughly 0-3 per cent. 


5. Details of the procedure in dealing with a nearly neutral system 
We suppose that the determinant C’ of the stiffness matrix has been 
aluated and found to be smail, so the system is near a state of neutral 

static stability. It is now convenient to consider an artificial state of the 

system which is exactly neutral and we can derive this merely by varying 
me element of C’: this element may be chosen so as to simplify the calcula- 
tion. The increment of the chosen element is C/K, where K is its co- 
factor. We next derive, merely by solving a set of linear equations, the 

4.10) neutral mode g, and may then proceed as explained in § 3. There is some 

idvantage in choosing the element of the modal column which is arbitrarily 

lumi kept fixed so that the product of this element and the varied element of C 


ecurs in the calculations 


6. Treatment of continuous systems 

The methods given above are not, in general, immediately applicable to 
rati ontinuous systems having infinitely many degrees of freedom. We may, 
L129 however, treat such systems with any desired degree of approximation by 
the ‘Lagrangian’ method in which the continuous system is represented 
by another having a finite number of degrees of freedom (2). This method. 
iso known as ‘the method of semi-rigid representation’, has been much 
ITO! used in dealing with aero-elastic problems and is of great utility. The 
TO! methods explained in the present paper can be applied at once to the semi- 
rigid representation of the system. 
rtl \nother method can be applied directly when the forces dependent on 
ne the displacements are derivable from a single-valued potential, for then 
The the stiffness matrix is symmetric. We saw in section 3 (6), that equation 


3.2) is valid to the first order of small quantities when (Cy is symmetric 
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whatever may be the nature of A and B. Now C, will be symmetric jn 
the circumstances we have postulated so (3.2) is applicable however many 
degrees of freedom the system may possess. To assist the application of 
the formula we now give a physical interpretation of the numerator and 
denominator. We know that the potential energy stored in the system 
when the displacement is qd, is 


V = 340 CQ. (6.1 
Next, the time rate of dissipation of energy in the system is 

D q Bq, (6.2) 
and if we suppose that q= qe" (6.3) 
this becomes D(A) = A*qy Bo. (6.4) 
Thus we see that qo Bdg = D1). (6.5 


where D(1) is the rate of dissipation of energy when 


q doe. (6.6) 
Accordingly equation (3.2) becomes 
oV 
‘ (6.7) 
D(1) 
The numerator and denorhinator can be calculated for a continuous system 
when the neutral mode of displacement is known. 
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ON THE STRESSES IN CLOSE-COILED 


n of . ‘ ’ 
? Pp ' 
HELICAL SPRINGS 
and 
ten By W. M. SHEPHERD (The University, Bristol) 
Re 17 Mareh 1950 
6.] 
SUMMARY 
The problen f the stress close-coiled helical spring is attacked by con- 
6.2 sidering the st1 es in an incomplete circular ring. It is found that the cross-section 
& 2 faring made from originally circular wire is almost exactly circular, but that othe: 
wo 
apes of cross-section would b ery considerably modified, especially if the ratio 
6.4) | f the coil radius to the wire radius were small. The stresses in a ring of cireulai 
6 ss-section are obtained and good agyeement is found with the results of Gohner. 
ov ° . ° 
ross-sections of non circular outline are considered with a view to reducing the 
iximum stress due to a given longitudinal force, but no worthwhile improvement 
66 the circular form is found to be possible. 
Introduction 
p Tue problem of the stresses in helical springs is one that has attracted the 
od . . . . 
nterest of engineers and mathematicians for many years. Many of the 
en ittempts to throw light on the subject have dealt with the problem of a 
close-coiled spring which, without serious loss 
| of accuracy, can be treated as an incomplete 
| : : > 
| ring or tore, subjected to loading as shown in Pp 
an Fig. 1. The first to make progress towards a so- \ 


ution was J. H. Michell (1) in 1900, who solved 
the problem for rings with oval cross-sections. \ | 


His method involved finding a function satisfy- 


V 

nga differential equation and taking constant | a, 
ae 
P 





alues on the boundary of the cross-section. An 





elementary and less rigorous solution was given 
by Rover (2) in 1913, and in 1931 Gohner (3), 


who appeared to be unaware of Michell’s work, produced a solution based 


Fie. 1. 


on the same differential equation. He, however, went farther than Michell 
ind obtained series solutions of the problem for a variety of cross-sections. 

In this paper the problem is attacked again, using the same differential 
equation and boundary conditions as Michell and Goéhner, but using a 
method of solution different from G6hner’s, being in a sense an extension 
of Michell’s method. The method is applicable to a wide variety of oval 
sections. 

Consideration is first given to the shape of the cross-sections of springs 
made from wire originally circular in section. An explanation is given of 


Yuart. Journ. Mech. and Applied Math., Vol. 3, Pt. 4 ( 
[Q J i Applied Vol. 3, Pt. 4 (1950)] 
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the surprising fact that the cross-sections of these springs are almost 
exactly circular. 

The spring of circular section is treated in detail and finally consideration 
is given to the possibility of reducing the maximum stress by using non- 
circular sections. It is shown that no worthwhile improvement is to be 
obtained in this way. 

Form of the cross-section 

Almost all helical springs in common use are made from wire of which 
the cross-section was originally circular. In the usual elementary treat 
ment of the problem it is assumed that, in the final form, the cross-section 
is still circular. This assumption is also made by those who have treated 
the problem more accurately. In view of the large permanent strains 
produced in the material during the manufacture of springs, it appeared 
necessary to examine this point carefully. 

Experiment shows that a straight bar of square section when bent to 
form a circular arc has finally a cross-section very far from square. It was, 
therefore, with some surprise that the cross-section of a valve spring taken 
from an aero-engine was found to be almost exactly circular. In this spring 
the ratio of the radius of the coil to the radius of the wire was about 8. 
A ground rod of the same section was then bent so that the corresponding 
ratio was about 3. Again the cross-section was found to be almost circular. 

The mathematical solution of the problem of the permanent strain in 
a bent rod involves considerable difficulty, but an approximate solution 
may be obtained which gives, in the author’s opinion, adequate explana 
tion of what was to him an unexpected result. 

When a rod is bent into a circular are the strain consists of two parts 
the elastic or recoverable part and the plastic or permanent part. On 
releasing the rod the elastic strains disappear almost completely, leaving 
only the plastic strains. The plastic part of the strain can be assumed to 
involve no volume change and will correspond to displacements of approxi- 
mately the type given by the equations of elastic bending if Poisson’s ratio 
is taken as 0-5. 

In Fig. 2 Ox, Oz are axes in the plane of a cross-section, Ox being in the 
plane of bending. Then the well-known theory gives the component dis 
placements wu, w. Apart from rigid body displacements u and w are given 


] . 
» u (z2?—a?)/4R, 


w 2 oh. 
where R is the radius of curvature of the strained centre line. These 
formulae are chosen so that the centre of the cross-section is not moved 


by the straining. It must, of course, be realized that these equations are 
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nly strictly valid if the strains are small. In the problem considered the 
strains may be by no means small and so the equations can only be 
expected to give an approximate picture of the form of the distortion. 


Fig. 3 shows the original and strained forms of the cross-section using 
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a ratio Ra = 4, a being the radius of the cross-section. In the same 
diagram are shown the strained and unstrained forms of the cross-section 
of a square bar deduced in the same way. Experiment gave very fair 
confirmation of these boundaries. It will be seen that despite the very 
considerable internal strains the boundary of the cross-section of the bent 
circular rod is, for all practical purposes, still circular. On this approximate 
theoretical basis the greatest deviation from the circular form should occur 
near the diameter at right angles to the plane of bending, but if R/a = 4, 
this only amounts to about 0-5 per cent. of the radius. Measurement of 
a bent rod indicates that the actual deviation is less than this theoretical 
value. For larger values of the ratio R/a the deviation will, of course, be 
smaller still. 

It appears that the circular section is the only one which preserves its 
form in bending. Even elliptical sections will be considerably distorted if 
the major axis is in the direction Oz and is large compared with the 
minor axis. 

This investigation casts considerable doubt on the usefulness of the 
investigations into the stresses in springs of square section. To obtain a 
spring of which the cross-sections are square, the form of the wire used 


would have to be carefully caleulated and would be by no means square. 


The twist of a circular ring sector 
Using cylindrical polar coordinates (7, 6, z), » denoting the distance from 
the axis of the ring, z the distance from its middle plane, and 6 the angular 
displacement about the axis of the ring, and following Michell (1) and 
Gohner (3), but modifying the notation, it is assumed that the only stresses 
are 6r and 6z and that these stresses are functions of r and z only. Then 
the only body-stress equation not identically satisfied is 
ero bz 2r0 
cr ’ Cz r 


and the compatibility equations give 


o7-0z 1cbz cOz 6: 


cr? Y or oz? y= 
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where ¢ is a function of 7 and z, and then equations (2) and (3) lead to 


ume 
blo} O-@ “h > ¢ ~ 
fai constant. ()) 
ull or= - r cr 
er : ° ° 
The stress 6n across the boundary is given by 
ent ‘ ‘ 
2-> Or- 
late Ay ré ~6. 
Cl Cs C8 
t vhere s. » denote the tangential and normal directions, and equation 
t of . 
then gives 
; 7 l od 
1a On > ° 
] 7™” OS 
Ie 
But Gn 0. so that 
t d constant along the boundary. (6) 
dif The problem then consists of finding a solution of equation (5) which 
the satisfies equation (6) on this boundary and which gives finite and con 
tinuous stresses within the boundary. 
the The method of attack used by Gohner involved finding a sequence of 
ina functions 4,, dy,..., each of which is constant on the boundary and such 
ised that the function 4, defined by 
are db 0; Ds PT he | d,,; 
pproaches a solution of equation (5). 
| The method used in this paper is in a sense the opposite of this. It 
ron a as ‘ ; é . 
onsists of finding a sequence of functions such that each is a solution of 
ulal “ as F , . . P 
' quation (5) and combining them linearly in such a way that the sum is 
and . . ? . 
onstant on the boundary. On putting 
ISISES | : 
h d b+ A+ B2*+ Cr 
" 
1 equation (5), 4, B, C being constants and % a function of r and z, it is 
| 
| found that —— Tr 
CW CW o cw 
0. 
or" Cz* rol 


Ifthe cross-section is symmetrical about Ox, then y% must be even in z and, 


if it is assumed homogeneous and of even degree 2n in r and z, then 


It follows without difficulty that 


2n—2p)(2n—2} Il jay, (2p-+2)(2p 2)ay,, 


nd that 7 As (), 


The functions of degrees 4, 6, 8 found in this way are 


; y4, 6r4z2?—7*, 16r4z4 


/ 


1676z2+-78. 
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It is easily shown that there are no functions of odd degree. Functions 
of higher even degree could easily be found if required, but for the purpose 
of dealing with the ring of circular section and with the other sections 
considered here, only the first two of these functions are necessary. 
The complete function ¢ is then taken in the form 
d A+ Br?+ C2z?2+ Drt+ 2F(6r4z?— 7) 7 


and the boundary is taken to be given by d = 0. 














x R-a 
Fic. 4. 
If the boundary of the cross-section crosses the axis z = 0 at points 
where 7 Ra, then equation (7) may be rewritten in the form 
PF {72?—( R—a)*}t( R--a)?— 77 re + hit KP + 62274, (8 


where A and k are two new constants which may be used to satisfy further 
conditions. For the purposes of this paper it is most convenient to pre- 
scribe the position and magnitude of the greatest width (in the z-direction 

of the section (see Fig. 4). 
For the circular section this greatest width must occur at + = R. 
La, so that regarding ¢ = 0 as the equation of the boundary, the 


conditions are 


... Ob 
(i) @ 0, (i) 0. when 7 R and z a. 
cr 


The second of these conditions gives 


2k = a*—18R?, (9 
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_— nd using this result the first condition gives 
ONE 
| 2h 52 R4— 20 Ra? +a‘. (10) 
Ons 
[he equation of the boundary now becomes 
: {72 (R a)y?(R t a)" rV18R2 a- 2r2) (11) 
" 1274+- 52 R4— 20 R2a?+-a4 : 
is now necessary to find out how nearly equation (11) represents a circle. 
SI | Ss iti iO! 
ipst ition Rx (12) 
| transfers the origin to the centre of the cross-section and the distance p 
| ofa point on the boundary from the centre is given by 
; 
p= xt 27, (13) 
On making the substitution (12) in equation (11) and using equation 
}), it is found, after some algebraic manipulation, that 
r2(a2— x?)(80R?+- 60 Ra+ 14a2—a?) 
ss p~ ‘ a 9.9 2 9 > ‘ 
6424+. 48 R3x+ 72 R22? — 20 R2a?+ 48 Ra + 1224+-a* 
i a of 6a? — 192? 
r?(a2—a?)| —_4 ' Lf. (14) 
4R? 16R! 
Since a—p is small the relation may be written 
pv-(a-— 27°) ‘ 
ad p > (15) 
8a R? 
proximately 
Equation (14) shows that the boundary fits the circle accurately at the 
i points whose 2 and z coordinates are (0, —a), (0, +a), (a, 0), (—a, 0). 
ints The greatest error will occur near the points where 2? sa*, and there 
p = 5a3/32.R? approximately. If R = 4a, the error is about 1 per cent. 
| for larger ratios of R a becomes rapidly smaller, being about } per cent. 
S } rr ° ° ° ° ° 
hen R = 8a. Thus it is seen that the boundary is circular for all practical 
t] purposes, 
then 
D1 Using the values of K and & given in equations (9) and (10), equation 
re 
Lins 8 becomes 
| 222)4 — 276 + 22(52 R4— 20 R2a?+-a*) —1r?(38 R4+ 30 R2a?)- 
i 


y4(22 R2+-3a?)+- 18 R6&— 37 Rta? 20 R2a4—a’. 


hen equations 4) give the stresses 


ee 22 
aad 2427? (52 R4— 20 R*a? +a’), (16) 
l - 2 
Q re 18227 +- 12) +(22 R?+-3a?)r += (38. R4+-30R2a?). (17) 


92,12 Hh 
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The resultant of these stresses acting on a cross-section must be equiva- 
lent to a force P acting along the axis of the ring. This force P is given 
by ae 
r | | 02 drdz 
taken over the cross-section. 

The boundary of the cross-section is so nearly circular that only a 
negligible error is introduced by taking the integral throughout the interior 
of the circle p = a. This point is referred to again later. The integrals 
involved are all elementary and the final result is 

P = wF\152R°+-44 R8a?— 15 Rat — (152 R44+-120a? R?), ( R?—a?)]. 
On expanding the square root in powers of a/R this gives 
P 64 "| 49 (“ ) 215 


Fra® a 128\R) * 1024 


(“ ) higher even powers of (a R)}. 
R 
(18) 

An estimate of the error introduced by integrating over the circle instead 
of over the exact boundary was made by evaluating the integral numeri- 
cally over the small region between the boundary and the circle in the 
case R = 4a. 

The error was found to be less than } per cent. and would be still less 
for larger values of R/a. 


The greatest stress occurs at the point r R—a, z = 0 and there 
A: q; r? = 0, where 
q , P P a 
; 128 R?a+ 160a2 R-- 1360 ( ; (19) 
k R 


From equations (18) and (19) the relation between the transverse force P 
and the greatest stress ¢ is obtained in the form 


; 5a ala *() e\~ 
qra2> 2R- ' 4R' 16\R 7 = 


r a ; 49 /a\2 215 /a\4 
“128\R} | 1024\R 


Gohner’s result, as given by Timoshenko (4), is, after an error in sign 





(20) 


has been corrected. 


F a\-! la l /a\? 

qra? on | r} '4R! a 

r a 7 3. (a/R) , 
16 1—(a/R)? 





The numerical values are set out in Table I. 
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from the circular form of the cross-section used by the author. This 
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TABLE | 


Table of Max. Stress « Cross-sectional Area/Transverse Force 


4 6 Ss 
hes I I 0 14°54 Io*7s 
| Snep! 10°95 14°79 15°09 

| ) 12 16 


For values of R/a greater than 8, the formula 
qra® 2R 5 
r ae '?2 
sof ample accuracy for practical purposes. 
The differences between the results obtained by Gohner and those 
btained by the author are not significant. They are due to the slightly 
lifferent method of approximating to the solution and to the deviation 


leviation is of the same type as the actual deviation, but is, as has been 


found, extremely small. 


Applications to sections other than circular 

The method used in this paper can be easily applied to the solution of the 
problem for oval sections other than circular ones. The shape of the oval 
section may be specified by the dimensions shown in Fig. 4. 

Taking the stress function in the form 

d fr? (R—a)*4( R+a)?—r?2 Hr? +k! 4+ K2?+ 622r4 

md the boundary to be given by ¢ = 0, the conditions at A and B are 
utomatically satisfied and it remains to satisfy the conditions at C 


nd D. i.e. 


d (). 0 when +r S. 2 b. 


these two conditions give equations to determine & and K. It will be 
oted that /& depends only on ¢ and not on b. 

The obvious next step is to find out whether it is possible to reduce the 
laximum stress corresponding to a given load by changing the shape of 
the cross-section. The common-sense method of reducing a maximum 
‘tress is to increase the amount of material in the highly stressed region. 
‘ince the greatest stress occurs at A it appears at first sight likely that 
section of the form shown in Fig. 4 with a b would be more efficient 
than the circular section of the same overall dimensions. In order to test 
this supposition, the problem has been solved for the series of dimensions 


R b C 
4 - 3°6, 3°/, 3°8, 3°9, 4:0, 4-1, 4:2, 4:3, 4-4, 
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the total load in each case being obtained by numerical integration. The 
results are shown graphically in Fig. 5, and it is found that the most 
efficient section is that in which the greatest width occurs at about 
} 4-2a, i.e. rather nearer to the outside of the ring than to the inside. 


Ss 





Pa 
we 


| 





Load 
Maz. stress X Cross-sectional area 























Ww 
aD 


38 40 42 44 


This unexpected result may perhaps be explained by the fact that the 
effective mean radius of the ring is increased or decreased according as 
the value of c is increased or decreased, and that the effect of this is 
sreater than the effect of the change of the distribution of the material. 
The greatest change in maximum stress possible is, however, so small that 
it is not worth while using a non-circular section for this purpose. 

The method put forward here could be used without difficulty to deal 


with the cross-section produced from an originally elliptical rod. 
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SOLATED FORCE PROBLEMS IN TWO-DIMENSIONAL 
ELASTICITY (11)t 


By M. ROTHMAN (Royal Technical College, Glasgow) 


Received 6 September 1949] 


SUMMARY 


In Part I of this paper the con:plex potential method due to Stevenson (1) was 
pplied to solve the problem of one or more isolated forces on the boundary of a hole 
nan infinite plate. In this part the same method is used to solve the problem of 

eral isolated forces in equilibrium on the boundary of discs of certain shapes. 


1. Isolated forces on the boundary of a circular disc with central 
fixing 

Ir there are a number of forces F, acting at the points A, (z = a,) on the 

boundary of a circular disc, then they must themselves be in equilibrium 





Fia. 1. 


r be balanced by external forces applied elsewhere. Suppose this is done 


by a force and couple nucleus /, Gp at the centre, so that 


K+>F=9, (1.1) 
Gt+hi > (F.a,—F,a,) = 0. (1.2) 


For the isolated forces on the boundary it is only necessary to assume the 


complex potentials: 


Q), -S Flogz,. z, = 2—4,, 
)- =" » ~~ ™ 
4 F log z, (F.d,— F_a,)log z,, 
TT —< TT — 
Forms part of a thesis approved for the degree of Ph.D. by the University of London. 


(Quart. Journ. Mech. and Applied Math., Vol. 3, Pt. 4 (1950) ] 
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and for the force and couple nucleus at the centre, Th 
2 F, log z 2xFj,zlogz , 2G, As 
Q, : Ws a) log z. 
‘ m7(«k+ 1) F a(k+1) rea so th 
The complete solution is then given by 
Q2 1) +2, 4+0,, w Wy tW1+Wy, 
' . A and ¢ 
where Q, and w, are expansible in positive powers of o, and Q and w 
satisfy the boundary condition 
Q(a?/z)+-(a?/z)Q'(z)+-w'(z) = k round 22 = a’. then 
_ - ~ ed 8. Zs 88 
Then as log Zz, = log(Z—4d,) lo : ) 
we have Di 
» ms 
(2, (2) + Q,(Z) + w4(z)+w4(z) - S F.‘log(—d,) + log z,— log 2} If 
TT homed r 
. a 
2K Ima 2— Fz . 
i ” _ (log a?—log z) 4 > Flog z,- a 
m(«+1) TT Ld To 2, 
—— ' a Kr 9, »; 
2x (44,—F,a,) 2h ee 21Gy and 
7 a Z m(k+1) 12 Also 
k’ 2G, ,2~: £4, 
1 T T . 
Tz 69 3 
9 
a* l 1\ _ 
Also | ja. ther 
22, iE. 2 
eo. , 2S ea 1 2a7F chec 
so that {OQ (2) +Q4(z)} 5 F,a,| ——, 
z 7 kam - dy m(«-+ 1)2° 
on :, a 2. ( 
hese equations reduce the boundary condition to T 
/_.9' 9 OW Al 
6 ; lx os iim 2a*h ; ths 
QQ, 1 —-Q)4(z)+-wo(z) 4 S (Fa,+F.a,)— ~: Pr’. than 
a eh ia, “aes sie (e+ | 22 
Z Z m2 n Z pen 
which can be satisfied by a suitable choice of k”, together with w, = 0 and 
and Q = Az+ Bz’, where A is real: for the boundary condition becomes \ 
Aa? Ba’ @ los oe 2a2k z 
|} ——_ +. _ (4 +2Bz)4 S (F.a,+ F.a,) ——~ 2 €', ther 
“a Z TT ham m(xk-+ 1)z* 
FR cire 
whence A - S (F.a,+ F.a,), B on. T 
2a? Ly a(x +1) 
the 
so that the complete solution is tha 
2 2F) log z l . 2,2? nis 
Q2 S Flog z, —— = S (F.a,+F,da,) + i (1.3) 
10 Lat a(«k-+1) 2rra*z Lon 7a*(«-+- 1) ma 
)- “ cm a ae >» nf ( 
22 ; 2 c-<— - 2«h,zlogz , 2iG 
Ww S F log z, { S (F.a, F’.a,)log i 0 i "log z. (1.4) ; 
- : - pare oe arlee © zer 
7 TT a m7(K-+-1) TT 
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ISOLATED FORCE PROBLEMS IN 


The pe riphe ral stress 


As before 2(nn+-ss) = Q'(z)4+-Q'(2), and nn = 0 round the boundary, 


30 that the peripheral stress ss is given by 


ann - Fz Fz) * (Fa, + F.a,) S - F : a 4 


7ma*(k+1) TA" ng T 2a? a,z—4a, 
LS = — Al 
. I f F(z—4,)+F(z—a,) 
F.z+F.z—(F.4,+ F.a,) 
then 
~ F\z+ Foz - Pa+Fa,) l > Fz- ate (F,4,4 Fay) 
ma*(k+1) 2a*n7 4y 7 ha 2a*—a,z—4d,2 
Diame trically opposite force S 
If there are two forces F, iR at z = ia, and F, = iR at z id. 
‘and G, being zero, then 
Fi a,+F, 4a, 2 Ra, F,a,+ yd, 2 Ra 
nd therefore S (Fia Fa.) (27a*) 2R (za). 
\iso = 9 
I Fi z—(F,a F, a,)}/(2a?—a, Z—d, z) Ra 
I F, F, d,+ F, a,)}/(2a?—a,Z—d,z) = Ria 
, 2R 2R 
therefore SS 0 0 


Ta Ta 


hecking with the well-known result. 


2, Curvilinear polygonal disc under isolated boundary forces 
The transformation z = co(1+Ao*"), where m is integral and greater 
han zero, and c, A are real positive constants, maps the interior of a 


urvilinear polygon on the interior of a unit circle in a complex o-plane. 


nd is everywhere conformal if z'(o) + 0 within the unit circle. 

Also as 2’(c) c|1+-(n+-1)Ao"| and z’(e) 0 gives oa” , 
A(n+ 1) 

hen for 0 n-+-1)A | there will be no critical points inside the unit 

ircle 

This transformation is of particular interest in two cases (i) when 7 2. 


the curvilinear polygon becoming a ‘dumb bell’ shape somewhat resembling 
that of the Standard British and French Cement Briquettes, and (ii) when 

is large and the disc also has central fixing, as we now have an approxi- 
mation to the problem of a cog-wheel (with shallow teeth). 


Case (i) In this case z co(1+-Ao*) and A }. Varying A from } to 


ero results in a succession of shapes ranging from a two-cusped epicvcloid 
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to a circle (Fig. 2), the ratio of the major and minor axes varying between 
2 and 1 in the extreme cases. (The cement briquette, however, though 
much of the same shape, has a ratio of 3:1.) 








Fic. 2. 


; = Re’? at the points z 
corresponding to o = o, = e'*r, then from the equilibrium of the external 


If there are a number of complex forces 


forces on the disc 
2,4, = 6, 2 4,% = 2 £4, (2.1) 


and the complex potentials giving these isolated forces are 


Q, a S F.log¢, where ¢, = o—o,, 
TT — 
a. AS cg > — _— 
w, S F log f+ ~ (F,.z,—F, z,)log ¢,. 
TT hmm TT — 


The complex potentials for the problem are then of the form 
Q) OQ,+Q,, WwW WoT 4; 
where 2,(¢) and w,(c) are expansible in positive powers of o and are 
chosen to satisfy the boundary condition 
2’ (0)Q,(0-1) +: 2(67)Q' (oc) +- a (0) +f(0) (0. 
where 


flo) = zo)! -= S Flog t= ae) 5 4 


, 


» mi os 2 gum ys Ly? - 
2x Fz 20 £2, fz, (2.2) 
i i 

As ¢,/¢, = —1/(,¢), then ¥ Flog @,/¢, is a constant; also 


(z—z,)/€, c(1+Ao?+-Ac,.a+-Ao*) 
and (Z—Z,)/¢, Z,0-1—Acao~*o-*— Aca; 3, 
whence (2.2) gives, after some reduction, 
}(c) (2c/7) > F (Ac? +-Ao,o) + (2¢/z7) > F(Ao-*0-?+-Aoo-3) + 
+(2/z0) > F.Z,, 


os 














ISOL 


so tl 


acco 


then 


as t! 


wh 
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twee 30 that f(¢) contains powers of o ranging from o*! to o-* only and we 
ough wecordingly take 

42, S a Wo c¥ Bo’, 
= T 
then 
2'(a)Q) (0-1) +-2(a0 7)Q,(c) +(e) 
cHh(A, A,) 3A(Az A,)|o 1 -(A,+ 2AA,)o~?-4 
(4,+AA,)o-3+ (3AA,4+ 24,4 B,)+(3A4,+-34,+2B,)o}; (2.3) 
is this must be equal to —/f(c) it gives five equations for the five constants 
1,, Ao, Az, By, By, namely, 
1,+4A,)+3A(A 1,) = (—2/ne) > FZ, 
1,+2AA, (—2A/z) > F.o7? 
1,+-AA, (—2A/z) > F.o-3 
ern as 
AA,+2A,+ B, (—2A/z) > F.o?, 
: 14,+34,+2B, = (—2A/z) > F.o,, 
whence solv Ing 
A ; S Fz 
1, - » Ua Fo?) a 
a(1— 3A") <4 mc(1 5A“) 
2A Y = 
{ , (2AF.0-?—Fa%) 
- (1 $A") 
| YA . = > | asa + 7 . = © 
A. SI a = [| A?(6A— 1) S (F.c, 3+ F. 03) 
7 hom Lar | 3dA*) al 
are A> (F.o7!+F.¢,)]. 
2A ’ ’ vac H2 
RB = VY AF.o72—(1+10A?) > Fo? |, 


Disc unde? fension 


If there are two equal, opposite, and horizontal forces 


F, T’ at oa, l, Fy T at op —l, 
then substituting in the above equations 
oT o—l1) 
QO, = A,o+A,07, w= B,c%, QQ, log i} 
T o+] 


27 (6X2—A—1) 2TXA—1) 


where and A, 
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so that the peripheral stress ss, given by 
88 real part of Q’(c)/2’(o), 
Z 47 A(1—3A) 
reduces to Ss . 7 
me (1—3A?)(1+ 9A?+- 6A cos 27) 
(see Fig. 5). This vanishes all round the dise for the two cases A 4 and 
A = 0, i.e. the two-cusped epicycloid and the circle (which agrees with the 
known result). 
Also at a point o = ia (a point on the mid-section of the disc) it is 
readily shown that 
vetyy = Q'(a)/2'(o)+-Q'(6)/2'(a) 
| A, +-34A,07—47'/n(o? 1)| [e( . 3Ao°) | 
2| —2T /74+-2Ta? +47 /7(1+a?)|/c(1—a?) 
(see Fig. 6). This is infinite at a +1 (the ends of the minor axis) and has 
the value 47'/zc at a 0 (the centre). 
Disc under briquette conditions 


If there are four symmetrically disposed isolated forces, as in Fig. 3, 


2 


F = Re-'P at o O71 em 
F, F, at o Os oa, 
i, = F, at oa 0}. 
and F, F, at o Gy e, *. 
A> RIR Ay 





ANR | RXAZ 





Then substituting in the above relations and simplifving, the peripheral! 
stress is found to be given by 


where 
38 4 Ri[(1—9A*){eos(a—B)+A cos(3a—f)}+ 12d? cos(3a+B)]|- 
+-cos 2y[ 3A(1— 3A)cos(a—B)+ 3A2(1 — 3A)cos(3a—8)- 
+- 18A2(1 —A)cos(3a+-8)]]/|ze(1 — 3A?)(1+- 9A? + 6A cos 2n) | 











and 


SS 


whic 
A 


mid. 


ther 
cor 


cou 


Th 


wh 


Als 
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and 
: 2 Ri sin(B+ y)sin(7+ «)+ sin(y—B)sin(y—«)| 
, mc(1+- 9A? +- 6A cos 2n)sin(» — «)sin(n + «) 
6A Ri sin(3y + «)sin(n—f)sin®(y-+ «)+sin(n+ 8)sin(3y—«)sin*(7—«) | - 
mc(1+-9A?+- 6A cos 27n)sin?(7— «)sin*(7+ a) 
which in the particular case A= 4, « = B = 4m further reduces to 


R xe, i.e. the peripheral stress is constant round the boundary. 
Also ea+yy = Q'(e)/2'(c)+0'(e)/2'(e), and at o = ia (a point on the 


nid-section) this reduces to (see Fig. 7) 


ve +yy t] — 18+-4(2a?—1)/(a*—a?+1)] | e(1—a?)]. 
Case (ii) In this case z co(1+Ao*") and A < 1/(n+1), and we suppose 
there to be isolated boundary forces F. Re! at the points Zz Z, 
corresponding to o o e'% in the o-plane), balanced by a force and 


couple nucleus at the origin given by F, and G,, so that 





K+ > F=0 Gotht > (F,z,—F,Z,) = 0. 
So A. J 
Yi -" 
) fo 
“ / 
\ 0" \ 
= 7 
: Y 
Fic. 4. 


The appropriate isolated forces F, at z, are given by the complex potentials 


a 
S2) = > Flog¢. where (, o—o,. (2.4) 
pS a 
w= —)> Flogt,+— > (F2,—F.z,)log{,, (2.5) 
-. TT — 
where Z c(o,+-Ao® +1) constant. 


Also the force and couple at the origin O are satisfied by the complex 


potentials , 
2F, log o 


: am(«+1) 


Z log Co l _ — = 
; S (F.zZ, -F.z,)log Go. 
‘ ) na? 
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Q = OQ, 


2’ (a){Q(a-1) 


Now, 


Q,(a0-4)-4 


and, 
Also, 
>... Bh) 
nF, 


a(Kk- 


(, 
I) 


w,(c) 1 w>(a) 


2z ( 


so that the boundary condit 


, 


z'(a)Q,(0-4) 
where 


1 sane 
(o)) —2 ¥ Flog 


a 
TT hmm 


f(c) 


2F, log o) 


m«+1){ 


TO Land 


sponding to w 


» 


T(c) 


TT << 


Now as z/o c(1+Ao”) and 
et Clo 

so that \2—2)/ 6, 

also (Z—Z,)/¢, 2/0 





+O, +Q,, w 


(c)log o4 


7 


LS (Fz 


+.) {FP (z—z,)—F(z- 
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= Wy TW, T 


The problem is therefore satisfied by the complex potentials 


Ws, 


where Q, and w, are to be expressed in positive powers of o only, and 
must be chosen to satisfy the boundary condition: 


k}+-2(a-})Q’(a) + wy’ (a) 0. 


») 


I 


)S Flogé, 
— 


ion becomes 


Z(a-!)Q(a)+- wy 


Sp k’ cae ~ > F. 
Co 


Td 


SS’ Plog f—k+ 


a} TO mmm 


y 2 F, logo 


m(k+1)’ 


(o)+-f(o) 0. 


loo ¢ 2«F, log a 
Og Cc. —+ oss 


a(k+1) 


(0 1)/ - F’ 2K, \ 2c Fi z 
| TT hand ¢, (k l)of zal(k+1)o 
F2)47 .F 2x £,2,—£,2, 
sities TT 1. ap Lat ‘ 


k’z, so that f(c) reduces to 


Zt Cp > 


- 
TOC —d 


Z/o 


The logarithmic terms cancel and the term in k’ can be neglected as it 
could always be eliminated by adding the rigid body displacement corre- 


(F.z,—F.z,)4 2(foxz— 2). 
m(xk-+-l)o 


c(o~2-+-Ao-"-?) then 
o,){1-+A(0"+1—o2+1)/(a—a,) 
z/o+-Ac(o,0"-1+ 02 o"-? +-o7") 
Ac(a; ly n+1_} o; 267” i oa; "a—*) 








[SOI 


and 
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and 
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ind so f(c) can be written 


2Ac 5 
a) \ Filia Co »o 2 om a; ) + 


2Ac ‘ . ‘ 
= % F (a> oe oC, “0 .. +o, "o-*)-4 


r 


» 
A 


Ss (eek: ~ LFye(1-+A0") + Fye(a-?-+Ao-"-*)], (2.8) 
y ) 


TT \k 


which the powers of o range from n to —(n+-2). We therefore take 


? 2 ” 

Q) S 1 o WwW c S B ot} 

= hos 0 hom r ? 

1 0 
nd the boundary condition becomes 
n---2 n+2 
c| | A(n l)o | S Ao c|a~ as * | > rA,o" 14 
i cy 
1 


c> (r+1)B,o"+f(e) = 0. 
1 


4 

t 
~ 
bo 


By +24 


tA(n+1)A,'-4 


ns 


NA, +A(n+1)(A,, 4 +444) } +o-7{A(n+- 1)A,,0+AnA, +A} + 


9° 


t 1 
oS (n+1)B (n+2)A,,,0} S {4,.+A(n+2—r)A, aa 
1 
> {(r+1)B,+(r+2)A,,.+A(n+1)A,_,}o" = 0. (2.9) 
1 


Comparing equations (2.8) and (2.9) we then have the following equations 
for the various coefficients 
2xAF, 


? T(K- 1)’ 
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and finally solving these equations, the A-coefficients are found to be 


A,[1—2A2(n+1)] ets (F.o.14+Fc,)- : > (F462), 
7 — TTC Md 
dW 9 ‘ ))\2 4 
A,(1—2nd?) ne (1 —Ain+-1))— FS Pe 2S Pont 
7 a(Kk-+- I a Tt it 
» y » m » y 2 | 
AnA, ants [1 A*(n+-1)] = S Fo 2x Fol] An 1) bs 
mK TT howd a(«+1)(1—2nA?) 
)»)\2 — » “ 
eS has a. 
1—2nd7| aw aw 1 Lod 
» ‘ é ))\3 i 4 
= = S Fo,4 2s (F24+F 2) DS (For1+Fa)), 
{ 2KAF, 
ii tics a(k+1) 
» ‘ 
A. = sa S {\(n+2—r)F.—F, ot +\o7" 
n| 1 A*(n +2 r)r| 


and substitution in the above equations will then give the B-coefficients 
if required. 

The cog-wheel 

If there is a single force +iF at the point o, 1, and a force —iF and 
couple G at the origin, we have the case of a cog-wheel with very shallow 
teeth under the action of a force and couple at the axle and a force at the 
actuated cog. 


In this case, 


— acm 
Q, ais 1), iz. 2iF log o, 
7 ig a7(K-+ 

21 Fz 41 Fce(1- 
Wy, : log(a—1)- = .; em 1), 
2xi Fz 21 Fe(1- 
Ws Kil loga Pe =] 1). 
~ a(k+ | 7 


Also substituting in the various equations for the constants: 


i 0. 4. 21 F[K(1 A)—A(I 7 is 2Ni 
= a(x+ 1)(1— 2nd?) 7 
j 2x F 4 2vi F[A(n+-2—r)+1] 
adi a(K- 1)’ ae 7 1—A*r(n+2 r)| ; 


i F[1—A2(n+1)] 2F . 2F[A(L—An)(e +1) +e(1—22n —A2)] 
a(k+1) z a(«+1)(1—2nA?) 


AnA 





~ 1 
“4T 


(xx + Jy) ( ne) 


) 


7v¢ 
Dp 


(x+y) 














and 
llow 


the 


(xcx+*+ YY 
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so that 


Q, = A,o+A,o7+ 
and the peripheral stress is given by 
— ova 
2.4 2.4 
real part of: | - 4. | 2A,0+nA, o"-!+4 
m(o—l1l)  xw(x+l1)o 


2NF 2xrF 
b= (n+1)o"+ —- (n 


2)or +1 


7  a(e+1) 


ro"-l(2NF Jar) “— 1+A(n+2—r) 
e| 1 L(n- 1)Ao” | Z 1—A?r(n+2—r) 


[ wish to express my further thanks to Mr. A. C. Stevenson for many 
helpful suggestions in connexion with the above problems. 
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